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CHAPTER i 


rtOt#5IW51 


t * §xJ 


OF POOR QUALITY 


INTRODUCTION 


The ablility tc predict the aeroelastic response oi 
aircraft wings is of increasing importance as attempts are 
made to reduce the weight of aircraft wings. One of the 
methods presently being explored to reduce weight is the use 
of a flutter suppression system (FSS) to reduce the required 
structural stiffness of the wing. The wing must be stiff 
enough to remain vibrational ly stable (positive damping) 
throughout ic-.- flight envelope. If the wing is not 
vibrationaV V ‘n «>ble. it will flutter and possibly cause the 


loss of the a!' 'raft. 

This pape' explains a method for predicting the 
vibrational stability of an aircraft with an analog active 
aileron F3S. Active aileron refers to the use of an active 
control svctea, connected to the aileron to damp vibrations. 
Uiing vit,' ^<t;ons are sensed by accel erometer s and the 
informatiD i used to deflect the aileron. Aerodynamic 
forces caused bi the aileron deflection appose wing 
vibrations and effectively add additional damping to the 
system. 

An assumed mode vibrational analysis approach is used 
with additional terms added to include the unsteady 
aerodynamics of a vibrating wing and the control system 
feedback. The assumed modes used are the actual vibration 
modes of the aircraft plus an aileron mode. The unsteady 
aerodynamic affects, modeled by the third order pade 
approximation method suggested by Edwards in the paper 
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"Applications o4 Laplace Trans-form Methods to Airfoil Motion 
and Stability Calculations", are used as the forcing 
functions for the vibration. 

A computer program called "SAFSS" was written to 
determine the vibrational stability of an aircraft wing 
using the method described above. The input information 
needed to use SAFSS" consists of; the natural frequencies 
of vibration, the generalized mass for each mode, the 
generalized force matrices for the mach number of interest, 
the control : aw matrices and the aileron deflections for 
each mode. ’SAFSS" produces a root locus plot from a Cal Comp 
plotter and a listing of the frecuencies and damping. 

A comparison between predicted and flight test data for 
the DAST ARW— 1 vehicle is made. DAST stands for drones for 
aerodynamic and structural testing. ARW--1 stands for 
aeroelastic research wing number one. The DAST ARW- 1 is a 
modified Firebee !I target drone fuselage mated to the ARW-1 
wing (see figure 1). ARW-! i? -- supercritical. sweptbacl. , 
transport-type wing with an aspect ratio of 6.B and a 
performance design point of mach 0.98 at 45,000 feet. The 
ARW-i wing is designed to be susceptible to flutter and is 
equipped with an active aileron flutter suppression system 


(FSS) . 


CHAPTER 2 


VlbRATION MODEL OF WING IN STILL AIR 

Sinco tliiitcr it a vibration problem, it would seem to 
be a reasonable idea to start with a basic vibrational approach 
and work up tc. the harder aspects one at a time. This section 
will review tne basic equations of the "assumed modes" method 
(sometimes called Rayleigh-Ritz method) used to solved vibration 
problems . 

An wing can be treated as a normal beam if no air 

is flowiTig over it. Because the aircraft i.s symmetrical, two 
types of vibration are possible; the symmetric case where both 
wings vibrate 180 degrees out of phase. The syimnetric and 
asymmetric vibrations do not couple with eachother so two separate 
analyses will oe run. The equations used to determine the 
frequem ies and damping of the system arc identical so no dist- 
inction IS made between the symmetric and asymmetric cases. 

Using the orthogonal normal modes (vibration modes) of the 
wing and generalized terms, each vibration mode can be written 
as a function of one variable and it's derivatives. The equations 
of motion are represented by 

( m . 1 •• u . ’ * + ( b . ) * u . * + ( k . ) ■“ u . = 0 
11 11 11 

Where "m^" is the generalized mass, "b^" is the generalized 
damping, "k^" is the generalized stiffness with "i" represent 
the i'th mode. Using the matrix notion: 






[M] = 0 0 . 

0 0 
0 0 • 

0 0 0 0 




the vibration equations can be written as one matrix equation 


[Bj U* = [K](ll^^= 0 

Lisinj; the fact that (U* is equal to I"), and reordering the above 
equation, the vibration equations can be written as: 



The left side of the abox'e equations can be simplified by 
multyplying by the inverse of the left hand square matrix. 
The inverse is 





-i 

Since CM3 is e diagonal matrix, CM3 is a diagonal matrix 
with terms equal to the inverse of the terms in CM3. 
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Multiplying both sides by the inverse solves the equations 
•for the derivatives of "U". 



Because both CM3 and CK3 are diagonal matrices, thei 
pi^ocjuct will also be a diagonal matri/%. 
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The same 


statement holds true -for the product o-f CK3 


and 


CB3. 
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Two relationships that are helpful in solving vibration 


problems are 


k/m =w**2 

i i i 


b /tn = 2*d #w 
i i i i 


where "w " is the natural frequency and "d^" is the damping 
factor. Lbstituting these rel ati onshi ps into the matrices 



a 




The natural frequency for a structure can be obtained from 
any program that does vibrational analysis such as "NASiRAN” or 
from a ground vibration test (a test in which the structure is 
shaken). Damping factors are harder to come by theoretically 
so a value is usually assumed. The wing pr ides only a small 
amount of damping so a small value, such as 0.005, can be assumed. 
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CHAPTER 3 


AERODYNAMIC INFLUENCE COEFFICIENTS 
The oscillatory motion of an aircraft wing in flipht will produce 
oscillatory forces on the wing. As the wing oscillates up and down, 
the deflection and its derivatives cause an effective angle of attack 
which changes the lifting forces on the wing. Likewise, as the wing 
oscillates torslonally, the changes in pitch and its derivatives will 
produce changes in t.he lift. 

The oscillation of the lift on the wing .lois as .t forcing function 
to the wing vibration and must be included in the vibration equations. 

The lift can he included in the followir.g manner 

IMl + IB] = fU 

whete '1.^ represents the oscillatory aerodynamic loads (general ired 
toms*. The aiirodynamic loads can be defined in the Laplace form to 
he 

L(s! = q " fQ(s)l * 

whore "q" is the dynamic pressure and [Q(s)l represents the aerodynamic 
influence coefficients. 

The [Q(s)j matrix is determined by curve fittisjg the aerodynamic 
influence coefficients (AICs) calculated at several vibrational frequencies- 
AICs for several reduced frequencies can be calculated by using a doublet 
lattice routine or some other unsteady aerodynamics routine. The 
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AICs arc a function of the reduced frequency and the mach number. The 
reduced frequency (k> is defined as 

k = (b * w; / V 

where "b" is the semi-span, ”w" is the frequency of oscillation and ”v" 
is the free stream velocity. A nondimensional form of [Q(s)l can be 
assumed to he 

-1 2 

lQ(s')J (sM - lku-1) * ([Pl'Ks') [P2’ls* + {F3M) 

where IPI* j. iP2' i. I P3" i and [RO* ] are constants deterrained by per- 
forming a least squares curve firt on the AICs calculated at the 
different reduced frequencies. The ”s’” is used to indicate a function 
of reduced frequency instead of oscillatory frequency. 

Kewriting this equ.'ition as a function of the oscillation frequency 
r e s u 1 ; s t : i 

-1 2 

igCsJj - tsr-b/v iKO' ) * (iPl • i(s*b/v) = [P2*3s--b/v + lp3'J) 

Multiplying [Q(s)J by U and labeling the product^X^ 

[Q(s)]\Ut 

The above equations can be rewritten by substitutins in for [Q(s)3 and 
rearranging 

(si + tRO*]*v/b) = ([ri'Ks)^ * b/v + [P2']s + [P3’3 * v/b>fu^ 

,'or the ease of vrriting and to agree with the nomenclature used by other 
authors the following notation will be used; 
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CR03 » v/b*CRO' 3 
CPi3 » b/vSCPi'3 
CP23 - CP2' 3 
iPZI = v/b*CP3'3 

Note that CPI 3, CP23, CP33 and CR03 must be recalculated for 
each different mach number. Using the above notion and 
taking the previous equation out of Laplace form, 

CI3CX'; - CP13CU’'> = -CR03CXJ + CP33fU3 CP23CU'> 

Remember inq that 


:m tU '3 + CB3 tU'3 + CK3 {U3 = <L3 


and 


?.L3 


q » CXJ 


the vibratJcn equations can be written as 
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Taking the inverse of the left hand side 
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and multiplying both sides by the inverse 


lO 


i 



r 



the vibration equations can be written as 



The above equation can be used to solve for the 
frequency and damping of the open loop aircraft wing 
vibration. 
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CHAPTER 4 


CONTROL LAW 

Although the procedure for putting e control law into 
matrix form is fairly standard, a few steps can be taken to 
simplify the solution of the vibration equations. In 
addition, part of the control law must be used elsewhere in 
the analysis (see chapter 5). 

A control law is defined as the output of a system 
divided by the input of a system. The input signal for this 
analysis is an accelerometer signal and the output is the 
deflection of the aileron. The numerator and denominator of 
the control law are generally written as the product of 
several first and second order Laplace polynomials (see 
figure 6). Part of the control law, a second order 
polynomial from the denominator and a constant from the 
numerator 

o 

g(s)=k/<s +c*s+k) 

where ”k" and "c" are constants, must be used for the 
aileron vibration mode (see chapter 5). The entire control 
law is equal to the product of its parts so 

G(s) = G’ (s) * g (s) 

where G(s> is the entire control law, G’ (s) is the control 
law without the aileron term and g (s) is the aileron term. 

The portion of the control law that needs to be put 
into matrix form is G'(s). Keeping in mind that the input to 



G*(s) is the accelerometer signal and the output is a state space vector 
the following equation can be written 

G* (s ' = V(s)/an(s) 

where "an” is the accelerometer signal and "V” is a state space vector. 
Taking G'(s) out of the Laplace form and putting it into state vector 
differential equation form results in the following equation 
^ ! A J Y + [G]^an3 
\V"^- !11J\Y^ 

where the matrices iAl. [Gj and [H] have no unique solution but depend 
on the state vectors 

In order to keep the vibration equations as simple as possible, it 

is desirable to have [G] and [H] contain as many zeros as possible. If 

the denominator of G'(s) is at least 3. orders larger than the numerator 

(Gj and [Ml tan !;o constructed to each contain only one non-zero tent. 

G‘(s) can I'c broken up into the product of its second order (or 

smaller) polynomials of the form 

2 2 

gHs) = (s + a*s +b) / (s + c*s +d) 

where g'(s) is a portion of G'(s), "a", "b”, "c” and "d" are constants 
n 

Assuming that g^*(s) is of the form 
B^’(s)=l/(s+d) 
and g^'(s) is of the form 
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2 

g^*(s) =1 / (s + c*s +d) 

2 

[G] and [H] can be forced to have a minimuni of terms. By assigning the 
state vectors^Y^to be the unlabeled inputs and outputs to each portion 
of the control law 

gj^'(s) = y^/ an 
= Yi/ Yj 

the following equations can be written 
y^ • =■ d*y j + an 
Y2’ ^ Yj 

y^' - Yj d*y,- cv^y^ 

Note that the accelerometer signal, ’’an" is used only once, meaning that 
[G] contains only one non-zero value. Assigning the state vectors in 
this manner also causes fH) to contain only one non-zero value. 





CHAPTER 5 


INCUil’ING CONTROL LAW IN VIBRATION EQUATIONS 
The movement ot theaileron will have two effects on the vibration 
equations. The first will be the addition of a mode shape (the aileron 
mode) and the calculation of AICs for the mode. The second will be the 
introduction of nondiagonal terms to IM] due to the addition of non- 
orthogonal aileron mode. 

In order to include the aileron as a vibration mode, terms for the 
generalized mass, generalized stiffness, and dampinf tern, must be used. 

In the previous chapter the control law was divided up into two parts 
u ( s ! - G ' ( s ) g ( H ) 

where G'(s) was put into state vector differential form and g(s) was of 
the form 

(’,< S / k ; ( S + b S + ki 1 

a * 

As staled 1 !; •!.< ‘act c;ia;.tc-i the output ol G’(s) is Since the output 

c rfs‘ IS ; Ji# .iiM-ion clef lection, the following equation must be true 
r'c \ 

where u is the aileron deflection. Substituting for g(s), to obtain 

k "VI-- u ’ ■ -t o "■ u ■ -t k * u 
ati^aaaa 
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Solving the above equation for ”u •*" 

& 

u * * = -k * u - b * u ' + k * V 
a a a a a a 

The vibration mode for the aileron can than be defined as having a 

generalized mass of 1.0, a generalized stiffness of k^, and a 

generalized damping term of b . 

a 

The input to the control law was described in the previous 
chapter as '.he signal from the accelerometer. The accelerometer 
signal is the acceleiution at the accelerometer location. Since 
the vibration model is based on superposition, the acceleration 
at any point, due to the k'th vibration mode, is equal to the 
acceleration state vector for the k'th vibration mode multiplied 
by the deflection k'tf mode shape at that point. Therefore, 



where "an" is the accelerometer signal, acceleration 

of the k'th vibrafioti mode and "r, " is the deflection at the 
accelorometci location lot the k’th mode shape. Ewpressing "an" 
in matrix mtatiori 

- lb I ^1.''^ 

where [D] is a row matrix containing the modal vibration deflections 
at the accelerometer location. 
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] 

A conversion factor may be needed between the units of the 
analysis and the units of the control law. If the control law is 
designed to convert a signal from g's (acceleration of gravity) 
to degrees deflection of the aileron, and the units being used in 
the analysis are inches and seconds, a conversion between the 
model and the real system must be made. Using the above example, 
the control law would have the units 

G(s) = AD f degrees) / an (gravities), 

h 

I and the analysis would require the units: 

2 

' G(s) - ua (modal deflection) c / an (in/sec ) 

where G(s) is the control law, "AD" is the aileron deflection, 

"u” is the state space notation for the aileron mode and "c" 

^ ,2 

IS the conversion factor. The conversion from g's to in/sec 

is simply 

1 g = (32.2 ft/sec^) * (12 in/ft) = 368.4 in/sec^ 
and assuming that a !0 degree •iefleciion of the aileron is 
equivalent to the aileron mode 

1 degree - !/10 deflect ion mode 

Therefore, the conversion factor is determined to be 
c = (1/368.4) * (1/10) 
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The equations related to the aileron mode and control law are 
\Y']= [A]^Y^+ [G]^an] 

iv3= fH]tY| 

‘^ai^= LDj1U“1 

u ••= -k * u - b * u * + k * c *\yl 
a aa aa a ^ 

Combining the above equations, the following state vector equations 

can be obtained; 

|Y'|= lA]^y^+ [G][D3tu*’^ 

U” ^-k*u-b*u’+k*c* [H]tY| 
a aaaa a ^ . 

The aileron mode is not orthogonal to the other vibration inodes so 
nondiagonal terms will be introduced into [HJ and/or [K]. In the case 
being examined here, there is no coupling between modes in [Kj but therc- 
is coupling in [M]. The terms in [M] are defined as 

where is the ele't.ent ij of the generalised mass matrix 

*’rj" are deflections at the point (x,y) due to the i'th anri j'th mode' 

shapes, and is the density of the wing at point (x.y;. Because 

the deflection of the aileron mode is zeio everywhere except the aileron- 

the integral will only be non-zere over the aileron. Assuming the aileron 

is rectangular and has a constant mass distribution, the generalized mass 

terms due to the aileron mode are 


- 


r^ (x,y) * r^ (x,y) ^-fCx.y) dx dy 
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m = -m * r * (r /3 + (r - r )/2) 
ia at il il it 

Where "ro " is the generalized mass of i’th row and the column rep- 
ia 

resenting the aileron mode; "m" is the mass of the aileron; "r " is the 

at 

deflection of the trailing edge of the aileron for the aileron mode shape 

the deflection of leading edge of the aileron for the i’th mode 

shape; and ”r " is the deflection of the trailing edge of the aileron 
it 

for the i'th mode shape. Assuming that the aileron mode is placed last, 
[M] would then be written as 


[M] = 


m 0 0 

1 

0 in 0 

2 

0 0m 


0 m 


la 


0 m 


2al 


0 m 


2a 


0 0 0 . . m m 

n na 

000 . . 01 . 

Putting the problem into state vector differential notation 
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The inverse of the square matrix on the left is 
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This i ■= the -form D-f the vibration equations that i*Ji 
sed solve -for the -freqaencv and damping cf 
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CHAPTER 6 


ANALYSIS OF DAST ARW-1 

The compute! program ”SAFSS" (Stability Analysis of Flutter 
Suppression Systems) was written to implement the analysis approach 
presented in chapters 2-5 (see figures 2 and 3). The input information 
to "SAFSS” consists ot the generalized force matrices, the control law 
matrices, the general i/od mass, natural frequency, and the aileron 
deflections for t-ach mode. In order to obtain the input information, 
other computet piogr.tms were used. 

The finite element analysis portion of the computet propran "NASTRAN" 
and an input liK provided by NASA Langley were used to calculate the 
first ten natui.ii frequencies and mode shapes of the DAST ARK-1. The 
aileron mode stiape was then substituted for one of the rif.id hodv modes 
calculated b\ "NASTRAN". The mode shajies were then input to the tioubl*. ; 
lattice puttie;. t •V-VSTRAN" and used to calculate t!u ,.ci .',d\ u., tit t 
influettcc coot f t' j.-ut 'Aits) at fifteen leducod trequenctes tseu 
A.), T.hc padi a!.ur>..xttiuiie curve fitting routine "QUEFIT" used the Aits 
at the reduced trequencies to calculate the gener.tlized fore*' matrices 
(see figure 5). 

The control law for the third test flight of DAST ARW-1 was obtained 
iron; NASA Langley atuS is shown in figure 6. The portion of tlie control 

law used for the aileron mode was 

2 2 2 
g(s) = (1256. 5) / (s = 502.7s + (1256.6) ) 
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Therefore, the aileron mode had a natural frequency of 
1256.6 hertz and a damping factor of 0.2. The remaining 
terns of the control law were put into matrices by the 
program "CONTRL'*. ■ . 

The input file documentation needed to use the computer 
program "SAFSS" is presented in appendiK B and a listing of 
the program is xn appendix C. A sample input file for the 
DAST ARW-1 symmetric, .825 mach number, 15,000 feet altitude 
case xs presented xn appendix D. The output for the sample 
case iS contained xn appendix E. 




CHAPTER 7 


FLIGHT TEST DATA 

Data from ttit tt.iid test flight of the DAST ARW-1 test vehicle was 
collected by the NASA Dryden Flight Test Research Center. The DAST was 
launched from a B-52 aircraft and remotely piloted to an altitude of 
15,000 ft, Testtrig was pertormed at the following mach numbers; 0.70, 
0.75, 0.775, 0. io and 0,825. 

The wing wa:- viorationally excited at each mach number to determine 
the vibrational abiiity. Excitation of the wing was produced by 
oscillating th. aileton in a continous frequency sweep from lu to AO hert.'-.. 
Excitation sweep:, w-'r* performed wtth the flutter suppression systen. 

\FSS) on a!td/or ot i (iiepend ing on the predicted stability at the test 
point) for both symmciric and asymmetric cases. 

The KSS used lot the HAST ARW-1 is active aileron control whic!; 
operate;, in tin ie-.l.-winp manner. Electrical signals fron tout accel- 
erometc! (tw. ;<■.■. it. d in i iic fuselage and one in each wing) are setit ;■ 
a c ompi'ttsati.; . 'Is. • ,.!«pei..-.ator separates the symmetlic. .isynmielr i c , .'iii 1 
rigid body motions then signals the actuators to hydraulically move the 
ailerons. A tim» history of the accelerometer signals is used to cal- 
culate the frequencies of vibration and damping factors. 

Because of at. eiror in the implementation of the control law, the 
Fss was operated at one-half of the designed gain. This error caused 
test data for a gain factor 
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of one-half to be compared with predictions for a gain 
factor of one. As a result, the DAST ARW-1 with the FSS on 
entered a flutter region near mach 0.B25 at 15,000 ft. and 
was lost. 




CHAPTER 8 


RESULTS 

The natural fre^iuencies for the synaaetric case of DAST ARW-1 obtained 
from the program "NASTRAN” and from a ground vibrations test (GVT), are 
presented in Table 1. For the frequency range of interest (10-<40 herta), 
the predicted and actual frequencies agree to within 0.5 hertz. The 
first wing bending node is of special importance because the aircraft 
fluttered in that mode. "NASTRAN" predicted the first bending mode to 
have a natural frequency of 9.1 hertz. The actual value was 9.6 hertz. 

The mode shapes and natural frequency were used to calculate the 
generalized force matrices therefore a 5* error could be considered 
large . 

The data taken at mach 0.755 is believed to be unreliable because 
the wing showed uninitiated oscillations, possible due to atmospheric 
turbulence. The KSS (flutter suppression sy; tem) was left on for 
velocities of mach 0.755 and above. The KSS off data for mach .0.755 and 
above was calculated from FSS on information. Due to an error in the 
implementation of the control law, the FSS on condition was operating at 
one-half the desired gain (K=.5) for the test flight. 

Table 2 shows the relationship between the predicted and experimental 
vibrational frequencies with the FSS offand Table 3 shows the relationship 
with the FSS on. The average error in frequency for the FSS off condition 
is 52 while the error for the FSS on condition is 82. 
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The FSS increased the vibrational frequency from 13.91 hertz to 19.93 
hertz. The analysis predicted the effect of the FSS quite well as an 
increase from IA.54 to 21.58 hertz. 

A fjraph of root locus versus mach number for the 15,000 ft. altitude 
case for the first wing bending mode is shown in Figure 7. Three gain 
factors are plotted: full gain (K=l, the designed gain factor), half 

gain (K=.5, tne actual Fss on gain) and no gain (K=0, the FSS oft gain). 

The root locus plot is interpreted by noting that as the real term 
approaches the imaginary axis from the negative real side, the system 
becomes less stable. As the iuaginary axis is crossed from the negative 
real to the positive real, the system goes from stable (no flutter) to 
unstable (flutter). The experimental results are shown to be consistently 
less damped and of a lower frequency than the predicted values. Trends 
seem to be predicted well, however more data is needetito draw any conclusions 

The predictions presented here are not the outcome of a single 
analysis, but rely on finite element, unsteady aerodynamics and vibration.il 
analysis wiiich lend to reduce the accuracy. The small difference wliic!: 
resulted can be caused by the inaccuracy (5*) of the finite elcroont 
analysis alone. 
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CHAPTER 8 
CONCLUSIONS 

A method for analysis of an active aileron control 
flutter suppression system has been explained and compared 
with flight data. The method was shown to produce reasonable 
results but relies heavily on finite element and unsteady 
aerodynamic analysis. The ablility of the finite element 
routine to match the mode shapes and natural frequencies of 
the ground vibration test will have a large affect on the 
accuracy of the flutter analysis. 

Future work in this area will move toward digital 
instead of analog control systems. Digital systems will 
reduce space and weight requirements as well as make 
possible the use of dynamic control laws. 
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TABLE 1 

NATURAL FREQUENCIES OF VIBRATION MODES PREDICTED BY NASTRAN 
AND MEASURED DURING GROUND VIBRATION TEST 



TABLE 2 

VIBRATIONAL FREQUENCIES OF FIRST WING BENDING MODE 
PREDICTED BY SAFSS AND OBTAINED FROM FLIGHT TESTS 
FOR FSS OFF AT 15000 FT. 





TABLE 3 


VIBRATIONAL FREQUENCY OF FIRST WINS BENDING MODE 
PREDICTED BY SAFSS AND OBTAINED FROM FLIGHT TEST 















Subroutine. 


K 



Figure 3. Flow Diagram o-f Main Subroutine of "SAFSS 
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(73a) (2250) » (a^ 76 // 8 S (^ 95 > 3 )^) (a^ ■» I20a ( 306 )^) (a^ ■» 100 s * (7lf) 

(s + 2) ({3 + 295.3)^ (3 + 1500)^ (3^ + 240s + (3^f2)^) (s^ + lOOs + (58)^) (s^ ♦ lOOs + (H2)^) 

(s^ 1003 •» (168)^) (295.3)^ (1^36.6)^ 

( b ^ + 76.70s + (295.3)^') (n^ + 589.4s + (439.8)^) (a’ + 502.7s + (J256.6)^) 


Figure 6 . Control hau for DAST A/V-1 Flight Test 3 * 
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APPENDIX B 


This appenJix contains a description of the input file for the 
computer proerair. "SAFSS" (Stability Analysis of Flutter Suppression 
System), This program is designed to use the output of an aerodynamic 
influence coefficient curve fitting routine has outlined by Edwards in 
"Applications of Laplace Transform Methods to Airfoil Motion and Stability 
Calculations". 

The input information includes the aerodynamics, the control law, and 
infornalion about each mode, such generalized mass, frequency of 
oscillation in stil! air and physical information about the wing. The 
input is grouped into four blocks. The first block contains information 
about the analysis options and the size of the problem. The second has 
all the information about the control law. The third block contains all 
the aerodynamic matrices and mach number. The last block contains the 
wing vibration and mode shape information. 

The term "carG'' means that all the inlormation is contained on one 
card. The term "card sot" will refer to a group of cards that contain the 
stated information. Card sets usually contain matrices. Matrices will 
be read in by rows (first row, second row etc.). The last value in a 
row of a matrix will be the last value read from that card. As an example, 
,< io!t by ton ;-:ntri\ to be input with only eight values read from each card, 
will need two cards for each row. The first card for each row will contain 


eight values with the second card containing the 



last two. 

Variable Name Column # Format 

Description 

CARD BLOCK 1 *« 

CARD 1: 

TITLE<6) 1-46 6AB 

The title and/or description desired on the data output 
and the root locus plot. The time and date will be supplied 
by the computer . 

CARD 2: 

NCASE 1-5 15 

The number at ditterent altitudes that calculations are 
to be pert or med at. 

MODE t - ! 15 

The number ot vibration modes Deing used tor the 
calculations. Must include the aileron mode. 

MALERN 11-15 15 

The number o-f the vibration mode that contains the 
aileron mode. Usually the first or the last. 

NGQibll 16-20 15 

The type of arithmetic progression desired for the gain 
factor, ( 1=1 inear , 2=geometri c) . 

Example; linear (0. , . 25, . 5, . 75, 1 . 0, 1 . 25, 1 . 5) 
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NGAINS 


geoisetric (0. , .25, .5, 1.0, 2. 0,4. 0,8.0) 

21-25 15 

The number of ^ain factors to be analyzed including zero. In the 
example above, noth linear and geometric have 7 gain factors. 

GAINUP 26-35 F10.6 

The desired value of the first gain factor after the pole (0.0). 

Both the linear and geometric series in the above example have a value of 
.25. 

lEGVEC 36-40 15 

If the eigenvpt-tors are to be printed, use the number "I". If any 
other number is round the eigenvectors will nor be orinted. 

CLEN 41-50 FlO.f. 

The span of the wing. It is used in the doublet lattice routine as 
the reference length to calculate the reduced frequency . 

IPLT 51-55 15 

If a root li.c ir plot is liesired. st t equal tc> "I". If any other 
number is useil, no nlut will b. made. 

YUkX 56-f : no. 6 

The maximum frequency (rad/sec) that is of interest. This value is 
used in plotting only, and does not effect the calculations. 
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CARD SET 3: 

ALT(I) 1-80 8F10.6 

This card(s) contains the altitudes to be used in the calculations. 
The maximum number ofaltitudes that can fit on one card is 8. If more 
than 8 altitudes are desired additional cards must be used. 



CARD A: 


NCONTR 1-5 15 

The size of the control matrix. 

CARD SET 5; 

AC(I.J) 1-80 4E20.13 

This set of cards contains the “A" matrix of the control law. It 
is a square matrix with dimensions NxN, wherein the order of the control 
matrix (does not include the terms used as the aileron mode). 

CARD SET 6: 

GC(I.l) 1-80 AE20.13 

This set of cards contains the "G" matrix of the control law. It 
is a row matrix with IxN dimensions, where N 
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is the order of the control matrix. 

CARD SET 7: 

HC(l.J) 1-80 4E20.13 

This set of cards contains the "H” matrix of the control law. It 
is a column matrix with Nxl diroinsions, where N is the order of the control 
matrix. 

CARD BLOCK 3 1 

This block of cards contains the PI * , P2', P3* , and R0‘ matrices 
as described by Edwards (ref. 1, 2, 3 and A). This block of cards is 
designed to be used as a unit. There is no need to separate the infor- 
mation in this block at any time. 

CARD 8: 

liACH 1-10 no . 6 

The mac!’, iiumber at which the calculations of the unsteady aer< - 
dynamics were performed at are included here. The unsteady acrodytuJUiU 
data is only valic for one Mach number. 

CARD 9: -- 

ISYM 1-5 L5 

This states whether the analysis is for a symmetric or assymotric 
case, where ”T" is for the symmetric and "F" is for the asymmetric case. 
These parameters have no effect on the problem, but does label the type 
of problem being solved. 

CARD SET: 
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PI(K.L) 


1-80 


4E20.13 


This set of cards contain the PI* matrix. This matrix is square of 
dimension MxM, where M is the number of mode shapes. One card is needed 
for every 8 elements in a row and there are M rows. 

CARD SET 11; ^ 

P2(K,L) 1-80 AE20.13 

Same as card set 10 except that it is the P2‘ matrix 

CARD SET 12: 

P3(K,L) 1-80 AE20.13 

Same as care; set 10 except that it is the P3* matrix. 

CARD SET 13:-- : — ^ 

R0(K,L) 1-80 4E20.13 

Same as card set 10 except that it is the RO’ matrix. 

CARD BLOCK 4 *******;;-*******Vrtt*^e.^*********«r*******iWir*irtt****-**^ 

This card block contains the physical characteristics of the wing 
and should remain the same throughout the analysis. 

CARD SET 14; 

OMEGA(I) 1-80 8F10.6 

This card contains the vibrational frequency of each mode shape. 

It is important that the order of the mode shapes is the same as the 
order for the doublet lattice routine. The frequency 
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for the aileron mode is obtained from the 2*nd order term of the control 


law and was not input in the '“Control** program. 

CARD SET 15;- 

ZETA(I) 1-80 8F10.6 

This card contains the damping factor for each vibration mode. The 
modal damping term for the aileron mode is obtained from the 2*nd order 
term of the control law that was not input into the ““Control** program. 

The modal damping term for the other modes were hard to obtain there- 
fore, a small value of approximately .005, was assumed. 

CARD SET 16: 

GKASS(I) 1-80 8F10.6 

This card contains the generalized mass terms for each mode. The 
value for the generalized mass term o' the aileron has a value of 
1.0000. If the user does not set the aileron term to 1.0000, the program 
will do so. 

CARD SET 17:-- 

DEFLECT(I) 1-80 8F10.6 

The deflection of the wing at the accelerometer location for each 
mode shape. The deflection for the aileron mode should be zero. 

CARD SET 18: 





XMD 


1-10 


F10.6 


The mass of the aileron in consistent units. If the problem has 
been done in units of inches, seconds, and pounds (as is usually the Case) 
the mass should be of such units as to be consistent with the rest of the 
problem. The consistent units for the above example are the mass of the 
aileron in slugs divided by 12 in/ft. The program does not do this type 
of conversion so tliat other system of units can be used. 

CONVFT 31-40 F10.6 

This card contains the conversion factor for the units of the input 
and output of the control law used in this analysis. If the control law 
input is in g’.*= (gravities) and the analysis is in inches and seconds a 
conversion must oe made. Likewise, the output of the control law output 
may be in degrees and the mode shape deflection may be one inch. For 

2 

this example, the input portion of the conversion is (lg)/(32.2 ft/sec *12 
in/ft). The output portion would be (1 degree)/(the angle due to the one 
inch deflect ion y . The proper conversion is obtained by multiplying the 
two terms. 

CARD SET 19:— 

The information below is required for each mode shape. The aileron 
mode will have a hinge line deflection of zero. Each card contains the 
information for one mode. 
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PHI (I) 


1-10 


F10.6 


The average deflection at leading edge of the aileron. 

Slgiil 21-22 FiO.6 

The average deflection at the trailing edge of the 
aileron. 
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<H)£CK SAFSS 

PROGRAM MAIN ( INPUT , OUTPUT , TAPE 1 » INPUT , TAPE6«0UTPUT ) 

c 

t: PROGRAM SAFSS (STABILITY ANALYSIS OF FLUTTER SUPPRESSION SYSTEM) 

C 

C« *#■»♦•»»■»#•» «■»«■»■»*>#■»«»■»««■■»•»»■«•*■»«•»»»*# 4* #*•»•*#•»«#*# #*#»»•»»» 

c 

C THIS PROGRAM PUTS TOGETHER AND SOLVES THE MATRIX REPRESENTING THE 
C WING MOTION AND THE CONTOL LAW. 

C 

C THE PROGRAM IS SETUP TO TARE THE CONTROL LAW MATRICES, THE 
C GENERALIZED FORCE MATRICES FOR A GIVEN MACH », 

a THE MODAL FREQUENCIES, THE GENERALIZED MASS AT EACH MODE, THE 

C: DEFLECTION AT THE ACCELOROMETER OF EACH MODE, THE AVERAGE DEFLECT- 

C ION OF THE AILERON AT THE LEADING AND TRAILING EDGE FOR EACH MODE, 

C AND THE MASS OF THE AILERON AND OUTPUT A ROOT LOCUS 

C PLOT OF THE FREQUENCIES OF VIBRATION FOR THE ALTITUDES OF INTEREST. 

C 

»«■»«««»■»*■)•*««««« *4* » 

C 

C DESCRIPTION OF INPUT DATA 

C 

C 

C TITLE = TITLE OF THE JOB 

C NCASE = NUMBER OF ALTITUDES TO BE EVALUATED 

C MODE - NUMBER OF VIBRATION MODE BEING USED 

C MALERN = WHICH VIBRATION MODE IS THE AILERON DEFLECTION MODE 

C lEGVEC DO YOU WONT THE lEGEN VECTORS PRINTED (0=N0) 

C NGAINF - NUMBER OF GAIN FACTORS TO BE SOLVED INCLUDING ZERO 

(.: NGAINT TYI't OF GAIN FACTOR PROGRESSION USED 

C (1 -LINEAR, Z'^GEOMEIRIC) 


C 

C 

C 


FIRST GAIN FACTOR TO BE USED AFTER ZERO 
ALTITUDE OF, THE I'TH ROUT LOlIUS PLOT 
SIZE OF THE CONTROL MAIRi;.' 



GAINUP 
ALT( I ) 
NCONTR 


n n n n n n 



- ; ■ 

f 

n 










ACn.J) « THE "A" MATRIX OF THE CONTROL LAW 

DCU.J) = THE "G" MATRIX OF THE CONTROL LAW 

HC<I,1» = THE "H" MATRIX OF THE CONTROL LAW 

P1(K,L) * THE S#»2 TERM OF THE FADE APPROXIMANT OF THE LOAD 

P2(K,L) = THE S TERM OF THE FADE APPROXIMANT OF THE LOAD 

P3(K,L) = -PHE CONSTANT TERM OF THE FADE APPROXIMANT OF THE LOAD 

RO<K,L) = THE E.,GEN'.<ALUE fERM OF THE PADE APPROXIMANT TO THE LOAD 

OMEGA! I )= FREQUENCIES OF THE CORRESPONDING VIBRATION MODE 

ZETA(I) = DAMPING OF THE CORRESPONDING VIBRATION MODE 

GMA3S = GENERALIZED MASS OF THE CORRESPONDING VIBRATION MODE 

DEFLECT = DEFLECTION AT THE ACCELEROMETER FOR EACH MODE 

XMD ^ NASS OF THE AILERON 

CONVFT = CONVERSION FACTOR BETWEEN CONTROL LAW UNITS AND THE 
UNITS USED WITH THIS PROGRAM. 

PHI(I) = AVERAGE DEFLECTION OF THE AILERON LEADING EDGE FOR 
THE i'TH MODE 

SIGd) = AVERAGE DEFLECTION OF THE AILERON TRAILING EDGE FOR 
THE I 'TH MODE 


DIMENSION W1 ( 1 1 , 1 1) ,W2( 11 , 1 1 ) ,W3U 1 , 1 1 ) 

DIMENSION PI ( 1 1 , 1 1) , P2( 11 . 1 1) , P3( 1 1 . 1 1) ,RO( 1 1 , 11 > , VP3( 11,11), 
1VR0( 1 1 , 1 1 ) ,DMEGA2{ 1 1 , 1 1 ) ,QBM( 1 1 , t 1 ) ,B( 11 , 1 1 ) ,VP1H2( 1 1 r 11 ) , 

2VP1 < 1 1 , 1 1 ) , T 1 TLE!8) , OMEGA! 1 1 ) ,ZETA! 1 1 ) ,GMASS! 1 1 ) , A ! 48 , 4B ) , WR ( ^8 ) 
3WI!48) ,Z! I 1 , 1 1 ) ,Z3!4a,48) , INT!48) ,ACC 15, IS) . AMTRX ! 33 . 33 ) , 
4BMTRX!33, 1 ) ,FMTRX! 1 , 1 ) , GMTRX ! 1 , 33 ) , HMTRX ! 1 , 33 ) , W4 ! 33 , 33 ) , 
f5PHI ! 11 ) ,SIG! 1 1 ) .REACT! 1 1 ) .DEFLECT! 11 ) ,ALT!20) , 

. GC ! 1 5 , 1 ) , HC ! 1 , 1 5 ) , W5 ! 1 , 1 1 ) , WB ! 1 5 , 1 1 ) , W7 ! 1 1 , 1 5 ) , Wa ! 1 5 , 1 1 ) , 

.W3!ll , 1 ) ,W10! 1 , 15) rWll ! 11 , 15) ,W12! 15, 15) 

'^•EAD aOl , !TITLE! I ) , 1 = 1 ,B) . 

CALL DATE!TITLE!7) ) 

CALL TIME!TITLE!8) ) 

F'RINT 901, TITLE 

READ 802, NCASE, MODE , MAlERN , NGA r NT , NGAI NS , GA I.^IUP , lEGVEC , CLEN , 

. IPLT.YMAX , 

READ 804, !ALT( I ) , 1=1 ,NCASE) 






READ B02. NCONTR 

MSIZE»M0DE#3+NCaNTR 

M0DE3*M0DE#3 

CALL PADJR (NCASE,«0DE.M0DE3.MSI2E,MAU£RN»CLEN.Wl.M2rM3. 

. P 1 . P2 . P3 . RO , y P3 . URO , OMEGA2 . QBM . B , UP 1 H2 » VPl , OMEGA » ZETA » GMASS , 

. A » WR , W 1 , 2 . 2*3 , 1 EOUEC , INT , AC , NCONTR , AMTRK , BMTRK , PMTRX , OMTRX , HMTRX . 
. W4» PHI ,Sir.,R£ACT, DEFLECT ,ALTrOC,HC.W5,W8»W7,W8/W9,WTO.Wl I #W12/ 
.NQAINT,NGAINS,GAINUP, IPLT, VHAK) 
aOl FORMA T(8A8) 

804 FORMAT ( OF lO.G) 

802 F0RMAT(5IJ5,F10.G. IS.FIO.B, I5.F10.6) 
aOl FORMAT! IHl . 10K,8A8) 

END 

SUBROUTINE PADJR(NCA3E*M0DF,M0D£3,MSI2£,MALERN,CLEN.M1 ,W2/ 

. W3 . P 1 , P2 , P3 , RO , i.'P3 , URO , 0MEGA2 / OBM , B , UP 1 K2 » UP 1 , OMEGA , ZETA f 
, GMASS, A, HR, WI, 2, 23, I EGUEC , I N T , AC , NCONTR , AMTRX , BMTRX , FMTRX , OMTrvX , 
.HMTRK,W4,PMI , S IG , REACT , DEFLECT , ALT , CiC , HC , W5 ,WG , W7 , W0 , WO , W1 0 , 

.HI I ,W12,NGAINT ,NOAINa,GAINUP, IPLT,yMA>() 

c 

(; GUBRQUTINE PADJR 

C 

#•»«■»«*#»####»# 

c 

C THIS PROGRAM RESI2ES THE ARRAYS TO FIT THE PROBLEM THEN READS 

C THE INPUT INFORMATION AND THEN GREATS AND 80LUES THE STATE SPACE 

C iMATRiX. 

C 

LOO I CAL IBYM 
REAL MACH 

DIMENSION PH I(MODE> ,SIG (MODE'/ , REACT (MODE) ,D£FL£CT(«ODE) / 

,WI (MODE, MODE) (MODE , MODE ), W5 (MODE , MODE ) , 

,PJ (H0D£',M0DE; ,P2( MODE, MODE) ,P3( MODE, MfJDE) ,R0( MODE, MODE) , 

. VP3 ( MODE , MODF; ) , URO ( MODE , MODE > , OMEGAZ < MODE , MODE > , GEM ( MODE , MODE ) , 

. a ( NODE , MODE i , UP 1 H2 { MODE , MDSE ), «.'P I ( MODE , MODE ) , 
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.aHEGA(MODE) , ZETA(MODE), 

.GMASS(MDDE) , A ( MS I2E , MS 1 2E ) , WR ( MS 1 2E ) , MI ( MS I2E ) » 2 ( MODE , MODE ) , 

. IMKMBIZE) ,Z3(MSI2E»M3IZE) ,ALT{MCAS£) /AC<WCOWTR/WCONT«) » 

. AMTRK(MQD£3/ MODES ) ,BMTRK(MODEa, 1 5 ,FMTf?M( i ,i) ,GMTf?Kn ,MUDE3) » 
.HMTf<X( 1 ,MDDE3) , W4 (MODES , MODES > ,GC(NCONTR/ J ) , HC (1 > NCOMTR ) , 

.WS( 1 ,MQD£> /WG(NCQNTR,MODE) , W7 ( MODE , NCONTR ) , W8 ( MCONTR y MODE > , 
.W9(M0D£y 1 > ,W10( 1 .NCONTR) .Wl 1 ( MODE , NCONTR ) , W1 2 ( NCONTR , WCONTR ) 
C»*nn>*»»*^**«*9***in****»*9******^m*ii-**'>i**i>«*i>«in>»«*»9*i-*»it**»*»*»inn>«» 

C READ MATRICES REPRESENT IMG CONTROL LAW (FROM PROGRAM "CONTROL") 
C*>'**‘>*<f^*»S‘0<^O*»'>****»1HHH>f»**(‘'^***9***^l>»0*»‘l>^'**1H>*****<lfHHHt»»****** 

DO 1 1 1*1 .NCONTR 

11 READ 801 , (ACC I ,J) ,J=1 .NCONTR) 

READ 801 . (GC( IR, 1 ). IR*1 .NCONTR) 

READ aOl . (MC( 1 . IC> . IC*1 .NCONTR) 

C READ MACH NUMBER. READ (T/F) WHETHER SYMMETRIC OR ASYMMETRIC 
C READ FADE APPROXIMATE MATRICES FOR UNSTEADY AERODYNAMICS. MATRICES 
n ARE PI '.P2'.P3'.R0', 

C4«ft«44«*4»«»»4««fl.(4«*»«4»«*fl*#«»4»il»<l«4»»44«4'll4««4*««»»4*»44###*#»4# 

READ 002, MACH 
READ (J03.ISVM 
1)0 12 IR=1,H0DE 

READ 00 1 . ( P 1 ( I R . I C ) . I C - I . MODE ) 

12 CONTINUE 

DO 13 IR* J .MODE 

READ 001 . (P2< IR. IC; . n> 1 .MODE ) 

13 CONTINUE 

DO 14 1R*I,MUDE 

BEAD B01,(P3nR.rC),IC*l.M0DCy 

14 CONTINUE 

DO IS IRsl.MODE 

READ 801, (ROCIR.IO.rC*!. MODE) 

13 CONTINUE 

C4444*4*-»4« 444444 »-;’<*»444«4»l(4»4.t44Ji»»»'!»««4444»4444«»4**44(!#4»4» 4 »»* 

C READ NATURAL FREQUENCIES OF NINE. OIBRATION. MODAL DAMPINO TERMS. 

C GENERALIZED MASS TERMS. AND THE DEflECTfON OF WINO AT THE 
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C ACCELEROMETER LOCATION. 

C»»##»*«**»»«*»*»-»*»*«*»***# ***«**»»*«*»»»»**»*«»*•»**<•*•** «»«**i»*«** 

READ a02. (OMEGA< I ) . I«1 .MODE) 

READ B02.(2ETA(I).I=l.M0Dt) 

READ B02.(GMASSm.I»l,M0DE) 

READ 802. (DEFLECTd ). I«1 .MODE) 

C READ T)I£ MASS OF THE AILERON. THE CONVERSION FAC < OR BETWEEN 

C THE UNITS OF THE CONTROL LAW AND THE UNITS OF THIS ANALYSIS. 

C*#*# »»»«»»»«♦««»«■*«**«««■»■«*'*»«»#«•«<»• »«»»*•*«»«»«■»***■*»■»»#**•«*•**»* 
READ B02.XMD.C0NVFT 

C READ THE MODAL AVERAGE DEFLECTION AT THE LEADING EDGE OF 
C THE AILERON AND THE AVERAGE MODAL DEFLECTION AT THE TRAILING EDGE 
C OF THE AILERON FOR EACH MODE. 

DO IG dl.MODE 
IB READ B02.PH1 ( I ) .SIG( I ) 

•*«««»*<««#■»*«* it 

C END OF INPUT REGION 

C PRINT INFORMATION FOR LISTING 

PRINT 904. MACH, MODE, NCQNTR.CLEN 
PiaNT 90G, <ALT ( i ) , d 1 .NCASE) 

PRINT 901 , (OMEGAC I ) , 1=1 .MODE) 

PRINT 902. <2ETAn ), 1 = 1 .MODE) 

PRINT 903, (OMASS( I ) , 1=1 .MODE) 

PRINT 900 
DO 17 IR=1,M0DE 

PRINT 907. < PI ( IR, IC) . IC= 1 .MODE) 

17 CONTINUE 
PRINT 909 
DO IB IR=1,M0DE 

PRINT 907, (P2< IR. IC) , IC- 1 .MODE) , 

10 CONTINUE 




I 


V/1 

CJ 



■> 


PRINT 310 
DO IS IR«l,f10DE 


PRINT 307, (?3( IR, IC) , IC=1 .MODE) 
13 CONTINUE 
PRINT 311 
DO 20 IR=i/mODE 


PRINT 307, (R0( IR, IC) , IC-^l .MODE ) 

20 CONTINUE 
PRINT 312 
»**»***»« 

C START PUTTING TOGETHER THE HATRIK FOR SOLUTION 


(J CALCULATE EFFECT OF THE AILERON MODE ON OTHER MODES 
C (NOND I AGONAL TERMS OF THE MASS MATRIX) 

SI=SIG(MALERN)-PHI (MALERN) 

DO 21 1=1 .MODE 

21 REACT! I ) = ->!M[>i«Sl4(PHI ( I )/2+(SIG( I ) -PHI < I ) ) /3) 
REACT(MALERN) =1 . 

GMASS(MALERW)=1 . 


*»#» INITA1,I2E PI..OTTER ** »»««*»« 
IF (IPLT.EQ.l) CALL INITPLT 


C START OF LOOP l OU EACH (iLIMUDE CALCULATIONS 
DO 500 !CASE=1 .NCASE 

C« 4*«4-K ■» 4 « *44 PLOT AC I S 4 4 4 * * * * 4 4 4 4 « 4 * 4 4 * 4 * 4 4 4 4 * 4 * 4 4 

IF (IPLT.EQ.l) CALL AX ISPLT < ALT ( NCASE ). MACH , YMAX . VMI N . 
.XMAX.XMIN.SCALE, IS YM, TITLE) 

C44444444444444444 4 *44*4444444 4 44*4*44444444444444444444*4444* 4 44444# 


C CALCULATE DYNAMIC PRESSURE AND AIR SPEED 

C4444«4*4«4*44444*»4444»»*44«»444*444S444444«44*44*«*#*444#*»#4#4**#4 

CALL QBARC ( MACH , ALT ( I CASE ) . U , QBAR . I NFLG ) 

PRINT 921 

PRINT 922, MACH, ALT ( ICASE), QBAR, U ' 

C 4 4444444*444*4 4444 *444*4444 4444 4# *444* *4*4 4 *4444 44 44444 »#44 444 4 #444 4 



n n n n n 



C: CALCULATE THE P3 AND RO MATRICES FORM THE P3 ' AND RO ' MATRICES BV 

C MULTIPLYINCi BY THE UELOCITY OVER THE SEMi-SPAN. CALCULATE THE PI 

C MATRIX BY DIVIDING THE PI ' MATRIX BV THE VELOCITY OVER THE 


’* C GEMI-SPAN. 

UB=U/<CLEN/2. ) 

CALL ARITH ( UB , P3 , 0 P3 , VP3 , MODE . MODE > 

CALL ARITH (UB.RO.O. ,RO,VRO,MODF .MODE) 

CALL ARITH ( UB , P 1 , 0 . . P 1 , VP I . MODE , MODE ) 

C CONSTRUCT THE K/M. B/M AND Q/M MATRICES. 

I DUm=OMEGA(MALERN) »»2 

D1IM2 = 2. IOMEGA (MALERN)«ZETA(MALERN) 

I DO 31 IR=1 .MODE 

I DO 30 IC«1 .MODE 

.! CtMEGA2( IR. IC) =^0, 

i B( IR. IC) --0. 

Vi QBM( IR, IC) =0. 

30 CONTINUE 

DME0A2(IR. IR)^0M£GA(IR)#*2 
G( IR, IR)=2.i»0MEGA( IR)»2ETA< IR) 

IF ( IR.NL.MAI I.RN) OBM( IR, IR) HIBAR/(iMAt;y( IR) • 

IF ( IR.NE.MALERN) QMEGA2( IR . HALERN ) «DUMH»REACT( IR)/GMASS( IR) 
IF (IR.NE.MALERN) B ( I R . MALERN ) ^DUM2<REACT ( IR ) /GMASS ( I R > 

31 CONTINUE 

construct: -P3 + pi«k/m 

-P2 + PUB/M 
RO - P1«Q/M 

CALL MULT < VP 1 .0MEGA2.WI .MODE) 

CALL ARl'THC 1 , ,VP3.~1 . .W1 .HI .MODE.MODL) 

CALL MULTC.'Pl.B.HZ.MCDE) 

CALL ARITH( 1 . , P2. -1 . .H2, WS.MOD) .MllDt ) 



CALL MULT(MP1.QBM,W3, MODE) 

CALL ARITH(-1 . ,VRD, I . ,W3»W3» MODE, MODE) 

XK=0. 

IDUM1=NGAINS+I 
DO 400 KK=1 , IDUMl 

C PLACE SMALL MATRICES INTO ONE LARGE MATRIX 

DO 32 IR=1,MSIZE 
DO 32 IC-1 rMSIZE 

32 A( IR, IC) >^0. 

DO 33 IR=I .MODE 

AIZ^IR-l + MQDE.ZttlR+MODE)^^! , 

DU 33 IC=l .MODE 

A ( 2» IR+MODE , 2* I C- 1 * MODE ) - -0MEGA2 ( I R , I C ) 
A(2«IR+MODE.2«IC+MODE)=-B( IR, IC) 

A( 2« IR+MODE. IC) =QBM<IR. IC) 

A( IR. IC) =W3( IR. IC) 

A( IR.2+IC-1+M0DE) -WI ( IR, IC) 

A( IR.2«IC + M0DE) =^W2( IR, IC) 

33 CONTINUE 

n-mnnnipinni it u mtmunniiimi^mnnniiiii »<i »*««««*•« 

C ADD ADDITIONAL TERMS TO TflE LARGE MATRIX THAT ARE CAUSED BV THE 
C ADDITION OF THE I LU7 I ElV SUPPRESS I ON SYSTEM. 

CALL FSS(A,0MEGA2.e,QBM,M0DE,MSIZE,NC0NTR.MALERN. 

, DEFLECT . REACT . GMAGS , VP I . AC , GC , MC , WTi . WG , W7 , W8 , WO . 

.WIO.WI 1 .W12.XK.C0NVFT) 

CfinnHnnnnnnniifitiiimii<iin>inmii-ii^iiii-ifii»-^iHHt<niiuniii)i>HHHnnnnntii»*»«it:ttnnnnnnni» 

C END OF MATRIX CUNSTUCTION. NOW GET THE EIGENVALUES 

CALL ELMHES ( MS 1 2E . MS I ZE . 1 . MS I ZE , A . I NT ) 

CALL ELT.RAN ! MS IZE , MS I ZE . 1 . MS I ZE . A . INI , Z3 ) 

CALL HORZIMSIZE.MSIZE. 1 .MGIZC. A, HR. WI .23. lERR.O) 

IF ZEROS WERE CALCULATED GO TO ZERO PLOTTER SECTION »*»*«» 
IF (XK ,EQ. 1000000. ) GO TO 42 





PRINT 931 ,>!K 
PRINT 932, lERR 

C#«*«*#* CALCULATE FREQUENCIES AND DAMPING **^^*********** 0 ********** 
DO 4 1 IE = 1 ,f1SIZE 
WNZ-WR( IE)*WR( IE)+WI(IE)*WI(rF) 

WN = SaRT(WNi!) 

2TA =-WR(IE)/WN 
C:yCLES=^WN» ,10034 

PRINT 933, WR(IE), WI(IE), WN . CYCLES, ZTA 
41 CONTINUE 

C PLOT RESULTS ON ROOT LOCUS PLOT. 

IF IIPLT.EQ.l) CALL POINT ( MS 1 2E , HR , W I , SCALE , KK , XK , YMAK , VMIN , 
.>:max,xmin) 

C««<nHHHHf* CHANGE THE GAIN FACTOR «»»««#««*«»«»««»«««•»«» »« 
IF (NGAINT.EQ. 1 ) XK=XK+GAINUP 
IF (NGAINT.EQ. 2) XK =GA INUP*2«’* ( KK- J ) 

IF (NGAINS.EQ.KK ) XK-1000000, 

GO TO 43 


C PLOT AND LIST THE ZEROS OF THE MATRIX 

42 PRINT 9S1 

PRINT 932, (WR{ I ) ,WI ( I ) , 1=1 ,MSt2C) 

ir (IPLT.EQ.l) CALL ZEROPLT ( MS I 2E , WR , N1 , SCALE , YM.^X , YM IN , 
XMAX,XMIN,NGAINS) 

C UNPACK THE E I GEN'.-'ECTOR MATRIX (2) PER EISPACK GUIDE 
C AND LIST. 

43 IF ( lEO'.'EC.CO.O) GO TO 400 
PRINT 941 

L.=-l 

in CONTINUE , 

IF (MSIZE.LE, 104L) GO TO 110 





L-L+1 
GO TO 1 1 1 
110 CONTINUE 

LL = MSr2:E-10#(L- l ) 

I. LL = 1 

113 CONTINUE ‘ 

I--S 

IF (LLL.EQ.L ) I =Ll. 

PRINT flGB 
I I = (LLL-1 )*1CM 1 

PRINT £142, ( (WR( J) ,WI ( J) ) , 1 i , H I ) 

F'RINT £168 

DO 51 J=1 .flSI^E 

ff?INT 942, (23( J.K) .K-n , I ir ) 

M CONTINUE 
UL^LLl-+l 

IF (LLI..GT.L) GO TO 1 15 
GO TO 1 13 
115 CONTINUE 

C ENDS LOOP. GO OACK AND CA1..CULATE NEXT GAIN FACTOR. 

•10O CONTINUE 

C ENDS LOOP, GO BACK FOR NEXT ALTITUDE 

a -It a-an-aa-rfsaata-aaattao .!»«»«« 

500 CONTINUE 
□ 01 FORMAT ( 4E20. i:i) 

802 FORMAT! OF lO.ij) 

□03 FORMAT (L5) 

004 FORMAT!/, 15X- "MACH - " , FQ . 3 , lO;: , I " VIBRATION MODES" 

. , //lOX, 12 , "TM ORDER CONTROL LAW" , lUX, "CMARACTERISTl C LENGTH-", 
,F7.3) 

.006 FORMAT!/IOX, "ALTITUDES TO BE E'.'AI.UATEI) , 2 ! / 1 OX ! I X , G 1 2 . G ) ) ) 
001 F0RMAT!/10>:, "MODAL FREOUENCIES' "./lO! IX.GIZ.G) ) 



302 F0RMAT</10X, "MODAL DAMPING " r / lOUX , G12 .S > > 

903 F0RMAT(/10>;, "GENERALIZED MASSES " , / 1 0 ( 1 X , G 1 2 . B ) ) 

908 EGRMAK /lOX, "PI MATRIX ••/) 

907 FGRMAT(10(lX,G12,e>) 

909 FORMAT! / lOX. "P2 MATRIX "/) 

910 FORMATI/lOX*, "P3 MATRIX "/) 

911 FDRMAT</10X, "RO MATRIX " / ) 

912 FORMAT! 1 OX." - END OF INPUT DATA _") 

9;-;i FORMAT! IHl ,5X. 10! " NEW CASE " ) / / 1 

922 FORMAT! lOX. "MACH = " , G 1 G . G/ 1 OX , " ALT - " , G 1 G . 13/ 1 OX . 

."GEAR =".G16..G/10X,"MTRUE=",G1G.G) 

931 FORMAT 1//20X, "THE GAIN FACTOR =".FG.3./) 

0;.i2 FQRMATl/lOX. "EIGENVALUES COMPUTED. ERROR CODE- " . 1 4 . / / 1 2X . 

. "REAL PART" , ax, " IMAGINARY PART " , 2 ( 9X , "FREGUENCV" , 3X ) , 

.4X, "DAMPING FACTOR", /yX.OlSX, "RAD/SEC". 5X> , ax. 
tn ."CYCLE/SEC"./) 

* 933 FORMA! !5X.4!F15.2.5X) .F15.4 i 

941 F0RMAT1///10X, "THE E 1 GENOEC'I 0R5 ARE",/) 

9G8 FORMAT!//) 

942 FORMAT! 10E13. 5) 

13S1 FORMAT! ///I OX, "THE ZEROS ARE LOCATED AT : " , / 1 2X , "REAL PART", 

.ax, "IMAGINARY PART" , /5X , 2 ! 8X . "RAD/SEC" , 5X ),/ ) 

9S2 F0RMAT!2!5X.Fi;>.2) ) 

RETURN 

e;nd 

SUBROUTINE CROSS ! A , B , C , I , J , K ) 

c 

C: SUBROUTINE CROSS 

C 

c 

c THIS PROGRAM MULTIPLES MATRIX "A" BY MATRIX "B”. MATRIX "A" IS 

C NECESSARILY AW I BY J MATRIX AND MATRIX "B" IS NECESSARILY AN 

C J BY K MATRIX, THE RESULT OF THE MULT I PL I CA,T I ON IS MATRIX "C" (I BY K) 

C * note: if the array "C" is the same ARRAY AS "A" OR "B" THE PROBLEM 



c: IS probable: screwed up 

c 

DIMENSION A( I . J) .B( J,K ) ,C< r .K ) 

DO 25 L=1 , r 
DO 25 N=1 
C(L»N)=0.0 
DO 25 M-1 . J 

C(L.N) =C(L,N)+A(L.M)*B(M,N) 

25 CONTINUE 
RETURN 
END 

SUBROUTINE INITPLT 

CiHHH*#*#*# *•»■»*■»**•« *»•»#*»##*«•»•»» «#!»»#■»•»■*»■»*#>«■»■»•» ■*###•)***«•»■»» 

C 

C SUBROUTINE INITPLT (INITIALIZE Pl.OTTER/ 

C 

C 

C THE ENTIRE 5 LINES OF THE PROGRAM IS SUPPOSED TO INITALIZE THE PLOTTER. 

C THESE INSTRUCTIONS ARE FOR A CALCOMP PLOTTER. IF VOU ARE USING 

C ANOTHER TYPE OF PLOTTER VOU MAY NEED TO CHANGE THESE INSTRUCTIONS OR 

C GET RID OF THEM ALL TOGETHER, THE OTHER PROGRAMS THAT DEAL WITH 

C PLOTTING INSlRUCnONS ARE A>! I SPI. I , PO i N f fi . AND ZEROPLT, 

C 

CALL PLGTS(0.0<4> 

CALL FACTOR! . 7074 > 

RETURN 

END 

SUBROUTINE AX I OPLT ( AL T , MACH , YMAX , YM I N t XMAX , XMIN . SCALE . I SYM . T I TLE ) 

C • 

C SUBROUTINE AXISPLT (AXIS PLOD 

C I 


c: 

C . THIS PROGRAM DRAWS THE AXISES AND THE IMFORMATIDN NECESSARY TO 
C DESCRIBE WHAT THE PLOT IS OP (MACH NUMBER » ALTITUDE. TITLE...). 

C OH, BY THE WAY THIS IS FOR A ROOT LOCUS PLOT. 

C THESE PROGRAM IS WRITEN WITH CALCOMP PLOTTER INSTRUCTIONS SO 

C TAKE IT OR VlRITE YOUR OWN. IF DO NEED TO CHANGE IT THE OTHER 

C SUBROUTINES THAT ALSO USE PLOTTER INSTRUCTIONS ARE TNITPLT, 

C POINT, ZEROPLT. 

C 

INTEGER TITLE(B) 

LOGICAL ISYM 
SCL^VMAX/11. 

NEXP=0 

IF (SCL.LT. 1 .) GO TO 30 
10 IF (SCL.LE.IO.) GO TO 50 
SUL = SCL / 10. 

NEXP=NEXP4l 
GO TO 10 

20 IF (SCL.GT..!) GO TO 50 
SCL-SCL»10. 

IMEXP = NEXP-1 
GO TO 30 

50 IF ( SCL . CiT . 5 . i SCALE 1 O .» 1 M . k <t NEXP 

IF (SCL. LE. 5. .AND. SCL.GT. 2- ) SCALE = 5 . -i 10 . »«NEXP 
IF (SCL.LE.2.) SCALE = 2.«10. '•"NEXP 

IF (SCL.GT.2.AMD.SCL.LE.2.5.AWD.NEXP.GE.!) SCALE-2. 5^« 10. «kNEXP , 
XMIN=“SCALE»G. 

XMAX = SCALEt»2 
VMIN^^O. 

CALL PL0T(0. ,0. ,3) 

C:ALL PLOT ( 10.735,0. , --3) 

IF (ICASE.EQ.l) CALL PLOT ( G . . 1 , . -3 ) 

CALL PL0T(0. , 1 1 . ,2) 

CALL AXIS(2. ,0. ,24H IMAGINARY AXIS ( RAD/SEC )i, -24 , 1 1 ., 30 ., 0 SCALE ) 
CALL AXIS(-G. ,0. , 19HREAL AXIS ( RAD/SEC ).- 13 , 0 ., 0 ., XMI N , SCALE > 



CALL SYMBOL(-G. ,-l .25. . 10, 12HGAIN FACTORS . 0 1 1) 

CALL SVMBOL<-G. .-1 .5, . lO.YHSYMBOLS.O. .7) 

CALL SYM80L(-B..12.,.14,11HALT!TUDE = 

CALL NUMBER 039. ,S99. .. H, ALT .0. .0) 

CALL SYMBOL 099. ,999. 14 , 12H MACH = ,0.,12) 

CALL NUMBER 029. ,999. ,. 14.MACH,0. ,3) 

If nSYM) CALL SYMB0LO99. ,999. . . 14, 14H SYMMETR-J C , 0 . , 1 4 ) 

IF (.NOT.ISYM) CALL SYMBOL 039 ., 999 .,. 14 , 15H ASYMMETRIC . 0 ., 1 5 ) 
CALL SYMBOL(-S, , > 1 .0, . lO.TITLE.O. ,90) 

RETURN 

END 

SUBROUTINE POINT (MSI2E, WR, WI ,SCALE,KK.>{K,VMAK,VMIN,XMAX,>{MIN) 

c 

C PROGRAM POINT (DRAWS POINTS ON THE PLOT) 

C 

C 

c; THIS PROGRAMS TAKES THE EIGENVALUES AND PLOTTES THEM IF THEM ARE 

C WITHIN THE AXIS OF THE PLOT. EACH GAIN FACTOR HAS A DIFFERENT 

C SYMBOL TO REPRESENT IT. "X" IS AWAYS USED FOR THE POLES AND "2" 

C IS SAVED FOR THE ZEROS WHICH ARE PLOTTED IN ANOTHER PROGRAM. 

C THE INSTRUCTIONS USED IN THIS PROGRAM ARE FOR A CALCOMP PLOTTER, 

C GO THEM MAY HELD TO BE CIIANOr’D. HIE OIIIER PROGRAMS THAT USE PLOTTER 

C INSTRUCTIONS ARE INITPLT, AXISPET- AND ZERGPLT. 

i:: 

DIMENSION WR(MGIZE) ,WI (MGIZE) 

IF (XK.EQ.O.) ISYMB=4 
DO 10 1:^1, MSI ZE 

IF (WR< I ) .LT.XMIN.OR.WRI I ) .ul .XMAX.OR. 

. Wn I ) .LT.YMIN.OR.WI < I) .Of .YMAX) GO TO 10 
XN = WR( I ) /SCAL.E 
YN=WI ( I ) /SCALE 

CALL SYMBOL ( XN , VN , . i 4 , I SYMB , 0 . , i ) , 

10 CONTINUE 



PX=-4 . 5+KK« . 6 

CALL NUMBER ( PX . - 1 . 25 , . 1 , >(K , 0 . , 2 J 
CALL SYMBOL < 999 . r 999 . , 25 , 0 . , - 1 > 

PX=PX+.3 

CALL SYMBOL ( Pa, -1 .45, ISYMB,0. , -t ) 

IF {ISYMB.E*Q.4) ISYM8--1 

I9YMB- ISYMB+1 

IF <IGYMB.EQ.4) ISYMB=5 

IF ( ISYMB.E(3.9) ISYMB = 9 

IVETURN 

L-ND 

SUBROUTINE 2ER0PLT < MS 1 2E , UR , W I , SCALE , YMAX , YM IN , XMAX , XM I N . KK ) 

C 

C: PROGRAM 2ERDPLT (PLOTTES ZEROS) 

• C 

C 

C THIS PROGRAMS TAKES THE EIGENVALUES NHICH REPRESENT THE ZEROS OF 

C THE MATRIX AND PLOTTES THEM IF THEY ARE WITHIN THE BONDS OF THE 

C AXISES. THE INSTRUCTIONS IN THIS PROGRAM ARE FOR A CALCOMP PLOTTER 

C AND MAY NEED TO BE CHANGED. THE OTHER ROUTINES THAT USE PLOTTER 

C; COMMANDS ARE INITPLT, AXI5PLT, AND POINT. THE ZEROS ARE REPRESENTED 

C BY A *'Z" ON nu: PLOT. 

C 

DIMENSION WR(MSIZE) ,WI(M3IZE) 

ISYMB=8 

DO 10 I-l.MSIZE 

IF (WR( I :• .LT.XMIN.OR.WRI I ) .GT.XMAX.OR. 

HI ( I ) .LT.'i'MIN. OR. HI ( n. GT. YMAX) GO TO 10 
XN'NRl I ) /SCALf: 

YN = NI ( n /SCAL.i' 

CALL SYMBOL ( XN, VN, .14, ISYMB, u. , - I ) 

10 CONTINUE t 

P:: = -4.5+<KK+1 .5)*^.G 


■< ( 7 ; 


f 

I 


CALL SV(180L(PX,-1 .20,24,0, . • 1 ) 

CALI. SYM80L< HX,-l .45,8,0. ,- I ) 

RETURN 

END 

SUBROUTINE FSS ( A , 0MEGA2 , 8 , QBM , MODE ,MS. ZE , NCONTR , MALERN , 
.DEFLECT,REACT,GMASS,OPl ,AC,0C,MC,W5,WS,W7,WB,W9. 

.W10,W11 ,W12,XK,CQNVFT) 

C 

C PROGRAM FSS (FLUTTER 3URPRESSI0N SYSTEM) 

C 

•»•»«» ft «•«#« if ftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftft 

c 

C THIS PROGRAM ADDS THE ADDITIONAL TERMS TO THE MATRIX THAT A FLUTTER 

C SURPRESSION SYSTEM CAUSES SUCH AS THE CONTROL LAW AND THE REACTION 

C OF THE MODES TO THE MOTION OF THE AILERON RESPONDING TO THE CONTROL 

C LAW. 

C 

C«*ftftftftftft*ftftft»ftftft»ifi«iHfftftftft#ft#ftftftiHfft4ftftftft»ft-i»lfftftft»ftft##*ft»*ftftftftft»*ft«ftftft*ftft-#»ftft*ft 
DIMENSION A(MSIZE,MSIZE) . 0MEGA2 ( MODE , MODE ), B ( MODE , MODE > , 

.OBM(MODE,MODE) ,DEFLECT(MODE) ,PEACT(MODE) ,GMASS(MQDE) . 

.OPl (MODE, MODE) , AC ( NCONTR , NCONTR ) ,GC( NCONTR, 1 ) , HC ( 1 , NCONTR ) , 

.W5( 1 .MODE) ,WS( NCONTR, MODE) , W 7 ( MODt , NCONTR ) , Wf; ( NCONTR , MODE ) , 

.W9(M0DE, 1 ) .W10< 1 .NCONTfT) , H 1 1 ( MODI: . NCUN I R ) . W12 (NCON IR , NCONTR ) 
N0DE3=M0DEft3 

C^ XK ft 0MEGA2 ( MALERN , MALERN ) «.JQN-.‘r r 
DO 11 L=1,M0DE 
W5( 1 ,L) -DEFLECT (L ) 

WS(L, 1 ) =-REACT(L ) /GMASS(L) 

11 CONTINUE 

DO 12. K=l, NCONTR 
N10( 1 ,K)^C»HC( 1 ,K) 

12 CONTINUE 

C»ftftft*ftft*ft CALCUl. ATE THE G^DftG.'M MATRIV! ftftftftftftftftftftftfttftftftftftftftftft 
CALL CROSS(GC, W5,WG, NCONTR, 1 .MODE) 

CALL CROSS (WG.GBM.Wa, NCONTR, MODE .MODE) 






DO 21 IR=l .NCONTR 

DO 2 I IC= 1 , MODE 

A( IR+M0DE3, IC) ^W0( IR, IC) 

21 CONTINUE 

CALCULATE THE -C.»D»K/(1 MATRIX #»*» a »<»»*#.»»*«.»« «« 
CALL CROSS (‘Wfi , 0NEGA2 , WS , NCQNTR , MODE . MOOE > 
no 22 IR*l,NCONTR 
no 22 IC=J,MODE 

A< lRff10DE3,2»IC-l-M10DE)-='-W0< IR. 1C) 

22 CONTINUE 

C« # a » « « a j» a « CALCUl.ATE TIiC -Ciaf)-*B/H MATRIU *«»»»•••»«» a » a <»■*»•*# # 

CALL CROSS ( WO, B, Wa, NCONTR, MODE -NODE) 

DO 23 IR“I,NCONTR 
DO 23 !C^1,MQDE 

A( IR + MODE3,2>*IC + MODE)«-wa< IR, !Ci 

23 CONTINUE 

»vM( CLACULATE THE C^Hl/M MATRIX » «# a « a » « » a*#* # a# » 

CALL CROGSIWO, W10,WI 1 ,MfJDE, I ,NCONTR) 

DO 24 IRal.MOm: 

I/O 24 IC=l,NCOHTR 

A(2»IR + M0DE, IC*-MaD£.l>----Wl U IN, IC> 

24 CONTINUE 

C»aa#a aaaaaaaaa CALCUI.ATE THE A a r,aI>-C«H/M MATRIX a a a ♦« * a #aa an 

CALL CRQSi)(WU,U I I , W I 2 . NCON1 R , MODE , WCON Tft ) 
no 25 iR=l,WC:iJNTR 
no 25 IC:=I,NCOWTR 

At IR-tM0Dt;:MCMlGDE3) -ACllft, IC) fW12< IR, IC) 

25 CONTINUE 

CaaiHf a a aaa a » a a CALCULATE THi; l-'laCaH/M MATRIX a a a ii a » a aa aa «* a a ■» 

CALL CROSSIMFI ,W1 1 ,W7,MODfc' , MODE , NCOWTR ) 

DO 26 I R= 1,1101/6 

nn 2 c, n>i,N(:cNiiR 

A1 IR, IC WIDDETJ) ^W7( IR, IC) 

2G CONTINUE 

RETURN , 

END 







/ 

' I 






SUBROUTINE ARITH (SA,A»SBrB.C,NR.NC) 

c 


C PROGRAI'1 ARITH < ARTIHMET 1C > 


C 

c 


c 
c 
c 
c 
c 
c 
c 
c 
. c 
c 
c 


THIS PROGRAM WILL MULTIPLE A MATRIX BV A CONSTANT AND ADD II TU SECOND 
MATRIX OF THE SAME SIZE AFTER THF SECOND MATRIX HAS BEEN MULIIPLED BV 


A SECOND 

CONSTANT. 




SA 

A 

CONSTANT THAT WILL 
NR BV NC MATRIX 

MHL 1 I PLV 

THE 

'■A'* MATRIX 

SB 

B 

CONSTANT THAT WILL 
NR 13 V NC MATRIX 

MULTIPLY 

THE 

"B" MATRIX 

C 

NR 

NC 

NR BV NC MATRIX (SOLUTION OF 
NUMBER OF ROWS 
NUMBER OF COLUMNS 

(SA<»A) + (GB»B)) 




DIMENSION A(NR.NC). B(NR.NC), C(NR,NC> 
DU 500 IRsl.NR 
DU 501 IC«1,NC 


C( IR, IC)^SA''Fi<. i(<, ICi t W , IC) 

•501 CONTINUE 
500 CONTINUE 
RETURN 
END 

SUBROUTINE QGARC <MACH» ALT.U.nOAR, INFLG) 

C» *•>»##«««!<»«»*>*<♦«» I*'** 1*, ■**»»«» 

c 

C PROGRAM QBARC 

C 

c , 

<; THIS PROGRAM lAKES THE MACH NUMBER AND AETI7UDE AND ESTIMATES THE 






C AIRSPEED IN FT/SEC AND THE DYNAMIC PRESSURE, 

t: 

REAL MACH 

DIMENSION ALTT(IO) , A(10>, RHOUO) 

DIMENSION ^MACH(10>,CLC0RRn0) 

DATA ALTT/0.0. 5000. # 10000. . 15000, .20000. ,25000, ,30CnM>, , 

.35000. ,40000. ,50000./ 

DATA A/1 1 16.45, 1097.09, 1077.40, 1057.35, IO'JG.92, 1016. H* 
.994.05,973.14,960.08,368.00/ 

DATA RHQ/. 0023769, .0020482, .001 7556, . 001 4962 ,. 00 1 2673 ,. 0010663 , 

. .00089069, .00073281 , .00058720, ,00036392/ 

D.TTA 0KMACH/.70, .725, .75, .775, .80, , 025 , . 05 , . 875 , . SO , . 925/ 

DATA CLCORR/1 ,025,1.026.1.027,1 .03,1.05,1.07,1.08,1.1,1.125,1. 155/ 

IF (ALT.6T.0.0) GO TO 102 

A1=A(1) 

RllOI =RHO( 1 > 

GO TO 103 
1.02 CONTINUE 

DO 200 1^1,10 
ISAM -I 

IF (ALT.LE.ALTn I) ) GO TO 10<< 

200 CONTINUE 
PRINT 930 

930 FORMAT <10X,**ALT IS OUTSIDE TARLtS. HILL USE RMfJ = 0.0, 01 = 968.07*) 

A 1=968.07 
RH0I=0.0 
GO TO 103 
100 CONTINUE 

D2=(Al.r-ALTT( IBAO-1 ) )/(ALTT< ISAO)-ALTT(ISAO-I) ) 

RHOI =RHO< ISAO- I )+02« (RHO< ISAM) -RHO( lGAO-1 ) ) 

AI “A( ISAO-I ) »D2*<A( ISAO)~A( ISAO-1 ) ) 

103 CONTINUE 
U=HACM^AI 

GBAR=u,5«RHai *U«U , 

C LIFT-CUROE SLOPE CORRECTION FACTOR 



n n 






IX=1 

DO 76 JA-1,10 

IF(MACH.GE.BKMACH( JA) ) IX=JA 
76 CONTINUE 

CORf-AC = CLCORR< IX) 
QBAR=GBAR»CbRFAC 
IF (INFLG.EO.O) RETURN 
tl^U*12 

Q8AR = Q£JAR/M4, 

RETURN 

END 

SUBROUTINE MULT ( A « B , C , I » 
DIMENSION A(I,I),B<I,I),C(I,I) 
DO 10 J-1 , I 
DO 10 K=1 , I 
XX^^O. 

DO 11 L-^1, 1 

11 XX-=>;X4A< J.L)«B(L,K) 

10 C(J.K)=XX 
RETURN 
END 


SUBROUTINE EL MtlES < NM » N . LOW , I CiH , A . 1 NT ) 423. 

C 424. 

INTEOER 1 , .1 . M , T4 . LA, NM , ) t.iU.h H 1 . L i)W , MMl , MF J 425. 

REAL A(NM,N) 

REAL X , V 

INTEGER INTdGHi 42S. 

430. 

THIS SUBROUTINE IS A TRANSLATION OF THE ALGOL PROCEDURE ELMHES. 431. 

C; NUM. MATH. 12, 349-360 < 1 9(36 » BV MARTIN AND WILKINSON. 432. 

C APPLIED MATHEUIATICS DIOISION, ARGONNE NATIONAL LABORATORV 465. 

C HANDBOOK FOR AUTO. COMP., VUL . I I -L I NEAR ALGEBRA, 339-358(1971). 433. 

C 434. 

C GIOEN A REAL GENCRAt. MATRIX, THIS SUBROUTINE 435. 

C REDUCES A GUBMATRIX SITUATED IN ROWS AND COLUMNS 436. 

C LOW THROUGH I OH TO UPPER HEHSENBERG FORM BY 437. 






C STABILIZED ELEMENTARY SIMILAR I TV TRANSFORMATIONS. 

C 

c ON input: 

c 

C NM MUST BE SET TO THE ROW DIMENSION OF TWO-DIMENSIONAL 

C ARRAY ‘parameters AS DECLARED IN THE CALLING PROGRAM 

C DIMENSION statement; 

C 

C N IS THE ORDER OF THE MATRIX; 

C 

C LOW AND IGH ARE INTEGERS DETERMINED BY THE BALANCING 

C SUBROUTINE BALANC . , IF BALANC HAS NOT BEEN USED, 

c SET Low=ii igh=n; 

C 

C A CONTAINS THE INPUT MATRIX, 

c: 

c ON output: 

c 

C A contains the HESSENBERG matrix. THE MULTIPLIERS 

C WHICH WERE USED IN THE REDUCTION ARE STORED IN THE 

C REMAINING TRIANGLE UNDER THE HESSENBERG MATRIX! 

e; 

C INT CONTAINS INFORMATION ON THE ROWS AND COLUMNS 

C INTERCHANfw^D IN THE REDUCTION. 

C ONLY ELEMENTS LOW THROUGH IGH ARE USED'. 

C 

C QUESTIONS AND COMMENTS SHOULD PE DIRECTED TO B. S. GARBOW. 
C 

LA IGH - 1 
KPl “ LOW t ) 

IF (LA .LT,. KPl) GO Ul 20>' 

C 

DO ISO M - HP 1, LA 

MMl - M - I ’ 


430. 

433. 

440. 

441 . 

442. 

443. 

444. 

445. 
44G. 
4 47. 
448. 
443. 

450. 

451. 

452. 

453. 

454. 

455. 

456. 

457. 

458. 
453. 
460. 
461 . 

462. 

463. 
464 . 
466. 

•* m: 

469. 

470. 
471 . 
4 72. 
473. 
4 74 . 



! 






X - O.OiiO 
I ‘ M 

DO 100 J - H, IGH 

IF ( ABS{A( J,MM1 ) ) ,LE. ABG(X)) 
>: = ACJ.MMl) 

I - J 

100 CONTINUE 



C 

INT(M) = I 

IF- (I .EG. M) GO TO 130 
C 

C: INTERCHANGE RONS AND COl.UliNG OF A 

DO 110 J - MM I, N 



Y - A < I , J ) 

A ( I , J > - A ( M , J ) 
A ( M , J ) - Y 


no CONTINUE 

00 120 J - n IGH 


Y - A ( J , I ) 

A( J. I) ^ AU.M) 

A ( J . N ) ^ Y 
1?0 CONTINUE 

C 

C END INTERCHANGE 

130 IF (X .EG. O.OEO) GO TO 100 

MPl ^ M ♦ 1 

C 

DO 100 I ^ fiPl, IGH 
y A ( n MM 1 1 

IF (V .£0. O.OEO) GO TO J RA 

Y Y / y, 

A(I,MM1) - Y 

C 

DO MO J ^ M, N 




475. 

47B. 

477. 

478. 

GO TO 100 

430. 
4B1 . 

482. 

483. 

484. 

485. 


487. 

488. 
483. 
430. 
431 . 

432. 

433. 

434. 

435. 
438. 
437. 


433. 
500. 
SOI . 

502. 

503. 
304 , 
505. 
50B. 
507 . 
SOB. 



n n n n ~ r; n o n n n n 



< 



.. \ 



o 


C 

C 

c 

c 


c 


c 


c 

c 

c 

c 

c 

c; 


140 A(I,J) = A<I,J) - V » A(M,J) 

DO 150 J ^ 1, IGH 

ir>0 A(J.t1) - A(JfM) + V I* A(J,n 

leo CONTINUE* 

ISO CONTINUi: 

200 RETURN 


509. 

510. 
511 . 

512. 

513. 

514. 

515. 
5IG. 
517. 
5 1 B . 


END 

SUBROUTINE EL TRAN < NM , N , L OW , I GH . A , INT , 2 J 

INTEGER I . J.N.KL.MM.HP.NM. IGH.LOW.MPl 
REAL A(Nf1. IGH) ,Z(NM,N) 

INTEGER INTI IGH) 

THIS SUBROUTINE IS A TRANSLATION OF THE ALGOL PROCEDURE ELfITRANS, 
NUM. MATH. IS, 1 8 1 -204 < 1 970 ) BY PETERS AND WILKINSON. 

HANDBOOK FOR AUTO. COMP., OOL . I I -L I NEAfy ALGEBRA. 372-305(1971), 

THIS SUBROUTINE ACCUMULATES THE Si.’.imiZED ELEMENTARY 
SIMILARITY TRANSFOIVNATIONS USED IN THE fvEDUCTION OF A 
REAL GENERAL MATRIX TO UPPER HESSENBERG FORM BY ELMHES. 

ON input: 

NM MUST BE SET TO THE ROW DIMENSION OF TWO-DIMENSIONAL 
ARRAY PARAMETERS AS DECLARED IN THE CALLING PROGRAM 
DIMENSION STATEMENT: 

N IS THE ORDER OF THE MATRIX; 

LOW AND Hill ARE INTEGERS DETERMINED BY T^E BALANCING 
SUBROUTINE BALANC. IF BALANC HAS NOT BEEN USED, 


520, 

527. 

520. 

529. 


531 . 

532. 

533. 

534. 

535. 
53G. 

537. 

538. 

539. 

540. 
541 . 

542. 

543. 

544. 

545. 
54S. 

547. 

548. 

549. 

550. 



SET LCW^l, rOH=-N; 


A CONTAINS THE MULTIPLIEf?S WHICH WERE USED IN THE 
REDUCTION BY ELMHES IN ITS LOWER TRIANGLE 
BELOW THE SUBD I AGONAL? 

t 

INT CONTAINS INFORMATION ON THE ROWS AND COLUMNS' 
INTERCHANGED IN THE REDUCTION BY ELMHES. 

ONLY ELEMENTS LOW THROUGH IGH ARE USED. 

ON output: 


Z CONTAINS THE TRANSFORMATION MATRIX PRODUCED IN THE 
REDUCTION BY ELMHES. 

QUESTIONS AND COMMENTS SHOULD BE DIRECTED TO B. S. GARBOW 
APPLIED MATHEMATICS DIMISION, ARGONNE NATIONAL LABORATORY 


INITIALIZE Z TO IDENTITY MATRIX 
DO 80 I ^ 1 , N 

DO GO J - 1 . N 
GO ZUrJ) ^ O.OEO 

ZdrI) ^ I. CEO 
80 CONTINUE 

KL ^ IGH. - LOW - 1 
IF (KL .LT. 1) GO TO 200 
FOR MP=IGH-1 STEP -1 UNTIL LOWH DO -- 


551. 

552. 

553. 

554. 

555. 
55G. 

557. 

558. 

559. 
5B0 . 
561 . 
5G2. 
563. 
5S4. 

565. 

566. 

567. 

568. 

569. 
5?0. 


572. 

573. 

574. 

575. 

576. 

577. 

578. 

579. 

580. 
581 . 


DO 140 MM = I, KL 
MP -• IGH - MM 


583. 

584. 


n n 







1 - INT(MP) 

IF- (I .EG. HP) GO TO HO 
C 

DO 130 J - MP, IGH 
Z(MP, J) 2(1, J) 

Z(I,J)-O.OEO 
130 CONTINUE 
C 

2 (I, IIP) = l.OEO 
HO CONTINUE 

C 

5'00 RETURN 
C 

END 

SUBROUTINE HGR2 ( NM , N , LOW , IGH , H , WR , WI , 2 , lERR , INUM) 

r: 

INTEGER I , J.K ,L,M,N,EN, 1 1 , JJ,LL,MM,NA,NM,NN, 

X IGH, ITS,L0W,MP2,ENM2. lERR 

REAL H(NM,N) ,WR(N) ,WI(N) ,2<NM,N) 
ivc AL p , G , R , i.i , r , w , , v , ra , a a , o i , or , 22 . norm . machep 

INTEGER MI NO 
I. OGICAL NOTLAS 
COMPLEX 23 
COMPLEX CMPLX 

REAL T3 ( 2 ) 

EGUU'ALENCE (Z3,T3(1>) 

THIS SUBROUTINE IS A TRANSLATION OF THE ALGOL PROCEDURE HQR2, 

C NUM. MATH. IG, 101-204(1970) BY PETERS AND WILKINSON. 

C HANDBOOK FOR AUTO. COMP., OOL . 1 1 -L INEAR ALGEBRA, 372-395(1971). 


590. 
591 . 

592. 

593. 

594. 

595. 
59B. 

597. 

598. 

599. 
GOO . 
601 , 

603. 

7. 

e. 

9. 


13. 

14. 


18. 

19. 

20 . 
21 . 
'T*'? 

23 . 


THIS SUBROUTINE FINDS THE EIGENVALUES AND EIGENVECTORS 
OF A REAL UPPER HESSENBERG MATRIX BY THE OR METHOD. THE 
EIGENVECTORS OF A REAL GENERAL MATRIX CAN ALSO BE FOUND 
IF ELMHES AND ELTRAN OR ORTHES AND ORTRAN HAVE 
BEEN USED TO REDUCE THIS GENERAL MATRIX TO HESSENBERG FORM 
AND TO ACCOmULATE THE SIMILARITY TRANSFORMATIONS. 

ON input: 

NM MUST BE SET TO THE ROW DIMENSION OF TWO-DIMENSIONAL 
ARRAY PARAMETERS AS DECLARED IN THE CALLING PROGRAM 
DIMENSION statement; 

N IS THE ORDER OF THE MATRIX; 

LOW AND ICiH ARE INTEGERS DETERMINE!) BY THE BALANCING 
SUBROUTINE BALANC. )F OAl.ANC HAS NOT BEEN USED., 

SET Low^^i , igh=n: 

H CONTAINS THE UPPER HESSENBERO MATRIX; 

Z CONTAINS THE TRANSFORMATION MATRIX PRODUCED BY ELTRAN 
AFTER THE REDUCTION 8V ELMHES, OR BY ORTRAN AFTER THE 

REDUCTION BY ORTHES, IF PlfRFORMtD. IF THE EIGENVECTORS 
OF THE HESSENBE.RG MATRIX ARE DESIRED, Z MUST CONTAIN THE 
IDENTITY MATRIX, 

ON OUTPUT; 

H HAS BEEN DESTROYED; 


WR AND WI CONTAIN THE REAL AND IMAGINARY PARTS, 

RESPECTIVELY, OF THE EIGENVALUES. THE EIGENVALUES 
ARE UNORDERED EXCEPT THAT COMPLEX CONJUGATE PAIRS 
OF VALUES APPEAR CONSECUTIVELY WITH THE EIGENVALUE 
HAVING THE POSITIVE IMAGINARY PART FIRST. IF AN 


ERROR EXIT IS MADE r THE EIGENVALUES SHOULD BE CORRECT 
FOR INDICES IERR+1 . . . , .N; 

Z CONTAINS THE REAL AND IMAGINARY PARTS OF THE EIGENVECTORS. 

If- THE I-TH EIGENVALUE IS REAL. THE I-TH COLUMN OF Z 
CONTAINS ITS EIGENVECTOR. IF THE I-TH EIGENVALUE IS COMPLEX 
WITH POSITIVE IMAGINARY PART. THE I-TH AND (I+D-TH 
COLUMNS OF Z CONTAIN THE REAL AND IMAGINARY PARTS OF ITS 
EIGENVECTOR. THE EIGENVECTORS ARE UNNORMALIZED. IF AN 
ERROR EXIT IS MADE. NONE OF THE EIGENVECTORS HAS BEEN FOUND*. 

lERR IS SET TO 

ZERO FOR NORMAL RETURN. 

J IF THE J-TH EIGENVALUE HAS NOT BEEN 

DETERMINED AFTER 10 ITERATIONS. 

ARITHMETIC IS REAL EXCEPT FOR THE REPLACEMENT OF THE ALGOL 
PROCEDURE CD IV BY COMPLEX DIVISION. 

QUESTIONS AND COMMENTS SHOUl.D BE DIRECTED TO B. S. DARBOW . 

APPLIED MATHEMATICS DIVISION, ARGONNE NATIONAL LABORATORY 


MACHEP IS A MACHINE DEPENDENI PARAMETER SPECIFYING 84. 

THE RELATIVE PRECISION OF FLOATING POINT ARITHMETIC. 
MACHEP = 1G.0£0»»(~10) FOR LONG FORM ARITHMETIC 
ON S3G0 

DATA MACHEP/ 1 .E-3/ 

DO 5 K=- 1 . IGH 
WR(K)-0. 

WI(K)=0. 

CONTINUE 

lERR - 0 . 


90 . 


n n r: n n 



SEARCH FOR NEXT E IGENUALUES 
(30 IF (EN .LT. LOW) GO TO 340 
ITS - 0 
NA = EN - 1 
ENM2 * NA - 1 


101 . 

102 . 

103. 

104. 


LOOK FOR SINGLE SMALL SUB-DIAGONAL ELEMENT 
FOR L=EN STEP -I UNTIL LOW DO — 

70 DO 80 LL - LOW. EN 

L = EN + LOW - LL 
f IF <L .EG, LOW) GO TO 100 

1 IP ( ABS(H<L.L-1 ) ) .LE. MACHEP « < ABS(H(L~1 .L-1 ) ) 

' X + ABS(H(L,L) ) > ) GO TO 100 

80 CONTINUE 
C 

C FORM SHIFT 

I 100 X = H(EN,EN) 

j IF <L .EG. EN) GO TO 270 

j Y = H(NA.NA) 

! N = H(EN.NA) 4 H(NA.EN) 

IF <L .EG. NA) GO TO 2.80 

If- (ITS .EG. 30) GO TO 1000 

IF (ITS .NE. 10 .AND. ITS .NT. 20) GO TO 130 
C 

C FORM EXCEPTIONAL SHIFT , 

T = T + X 


107. 

108. 
108. 


1 12 , 


114. 

115. 

116. 

117. 

118 . 

119 . 

120 . 



/ 



122 


n n n n 








t- 123. 

DO 120 I = LOW, EN 124. 

120 H< 1 , n - H( I , I ) - X 125. 

C 12B. 

a = ABS(H<EN,NA) ) + ABS < H < NA , ENM2 ) ) 

X = 0.75E0 * S 12B. 

Y = X ‘ 129. 

W = -0.437SE0 » S S 130. 

130 ITS = ITS + 1 ■ 131. 



LOOK FOR TWO CONSECUTIUE SMALL 

BUB-DIAGONAL ELEMENTS. 

FOR M=EN-2 STEP -1 UNTIL L DO — 
DU 140 MM L, ENM2 
M = ENM2 *• L - MM 
2Z = H(M,M) 


P = (R - b - W) / H(MH,M) + H<M.M+1) 
Q = H(M+1 ,M+1 ) - 22 - R - S 
,R = H(M+2.M+1 ) 

S=-ABS(P)+ABS<Q)>ABS{R) 

F = P / S 
Q f G / S 


133. 


135. 

13B. 

137. 

133 . 

139. 

140. 

141. 

142. 


144. 

145. 


R - R / b 14S. 

IF (M .EG. L) GO TO 150 147. 

,IF { A6G(H(M.M-1 ) ) » ( ABS<G) + ABSIRU .LE. MACHEP * ABS(P> 14B. 
X * ( ABS(H(M-1 ,M~n ) + ABS(2i:) + ABS(H(M-+1 .M+l> > ) ) GO TO 150 149. 

140 CONTINUE ISO. 

150 MP2 M f 2 152. 

Df) IGO I MP2, EN 154. 

H(I.!-2) O-OEO 155. 

IF (I .EG. MP2) GOTO IGO 15S. 

H<I.I-3) O.OEO J57- 

160 CONTINUE 15B. 

C DOUBLE GR STEP INODLOING ROWS L TO EN AND, 






c COLUMNS M TO EN 


DO 260 K = M, NA 

NOTLAS « K .NE. NA 
IF (K .EG, M) GO TO 17 (‘> 

P = H(K,K- 1 > 

Q = H(K+‘l ,K -1 ) 

R = O.OEO 

IF <NOTLAS) (? = H(K + 2 .K- 1 ) 

>; - AbS(P) + ABS(Q) + APSCR) 

IF (X .EQ. O.OEO) GO TO 760 
P = P / X 
Q = C / ;( 

R = R / X 

® = SIGN( SQRT(P»P+Q»Q»-R#R) ,P) 

IF (K .EG. M) GO TO IBC 
H(K,K-n - -B # >; 

00 TO 190 

ISO IF (L .NE. M) H(K,K- 1 ) - -H(!.,K- 1 ) 
190 P - P + s 
X = P / S 
V = Q / B 
22 = R / S 
Q = Q / P 
R - IV / P 

C , 

C ROW MODIFICAriON 

DO 210 J K. N 

P“H(K,J)-tQ*H(K<-l.J) 



IF (.NOT. NOTLAS) GO 

TO 

2 CO 


P “ P + R « H(K+2,J) 
H(!< + 2.J) = H(K + 2,J) • 

• P 

» 2 .. 

200 

H(K + 1 , J) .t H(K+1 , J) 
H(K , J ) ^ IKK , J ) - P ^ 

• P 
' X 

# V 

210 

CONTINUE 




J ^ MtN0(EN,K+3) 





162. 

163. 

164. 

165. 

166. 
167. 

169. 
1 70 . 

171. 

172 . 

174. 

175. 

176 . 

177. 
176 . 

179. 

180. 
161 . 
1 . 12 . 
133. 


135. 

186 . 

187. 

188. 

189. 

190. 
. 91 . 
192. 
133. 
194. 




c 

C CQLimw HOD I r I CAT I nw 

DO 230 1 *1,,/ 

P * X * / ! ( I , K ) ♦ V H < I , K ♦ { ) 
IF (.HOT. MOTl.AS) GO TO 270 
P ~ P ^ ZZ # Hil.K'hZ) 

MM ,K *2/ * H( I ,K>2) P « h 

220 H(l,KH1 * P « A 

H(1 ,K / * H(I ,K ) " P 
730 continue: 

c 

C ACCUflULATE ( RAN3F0RMATI0NS 

00 250 I * LOW, IGM 

P - X ,* 2 <1 , K > + Y * Zn.KH) 
IF (.NOT. WOTLAS) GO TO 240 
<» F ^ F ZZ # Z( I,K + Z) 

2<I,K + 2) “ Z(I,K-«-2) - P tt R 

240 2(1, K+I) = Z(I,K+1) - P « a 

2(1, K) ■=• 2(1, K) “ P 
250 CONTINUE 

C 

200 CONTINUE 

C 

GO TO .70 

n 

0 ONE ROOT FOUND 

2/0 (!<£N,EN) ^ H f T 
PiR(F:N> . H(£N,£N> 

WKEN) > O.OFO 
FN l</A 
GU TO CO 

(; 

(; TWO ROOTS FOUND 

;?00 P (V - 4i / 2.0E0 
0 - P « P f w 
22 - GORT( AD3(GI)) 




ISB. 
107. 
13B. 
ISO. 
ZOO. 
20 1 . 
202 . 
203. 


205. 

20S. 

207. 

200 . 

209. 

210 . 
211 . 
212 . 

213. 

214. 
213. 
216. 


218. 
210 . 
220 . 
221 , 
222 . 


224. 

223. 



•BfiZ 

*( 3 G 2 


*CS2 
•2:52 
• TS2 
•0S2 
*GlrZ 


• Z , fr 2 

' St’S 

• efr 2 

* frfr 2 

•GfrZ 


• r ^2 

•OfrZ 

* 8 CZ 

'^.CZ 

' 9 E 2 

•SE2 

* 1 ^ E 2 

•CEZ 


•OEZ 

*G2Z 

*022 

•LZZ 



r 

j ! 


*■ <) - (nb'DZ * a ■ <n3'1)Z 
'HJ'DZ *> t> t 72 * n - iyW'DZ 

■’ 22 

mn 'Mrn « i otc cm 
sNotiwuiu.mmfn ~'i.vmurrj'j^ d 

D 

annuNOD ooc 

22 » ii ■ •• 0 ' (N3'1)M i 

iH3‘ I )H 0 6 + ZZ r D - ' VN ' I ) H 

(t/W'IJH - 22 
N3 'I -- I OOC 0(1 
N()iiv3i jiaoH Nwmoo 0 

:j 

jnmiNOD OHZ 

2Z •> << ' # 0 ' (( "N3)H ^ 

(f'N3)M « d ♦ ZZ « 0 ' (T'ornw 

(r'OW)H ^ 22 Vi ,' 

N 'On - f OGZ IK! 

W0UV3T3HJ0W «0-iJ 0 

n 

/ 22 « O 

« / X = ,i i 

(ZZeZZ-tKi^HnmS -r n I 

I,*JN‘H3)U ^ y. I 

030*0 ^ (N3)IM I 

030*0 3 {ON)!;'; ! 

22 / M " (W3)8M (030*0 *3N* 2Z) 31 I 

(yN)UM (N3)(JM 
22 + X - (VN)iJM 
(d'ZZ)MOJfi d 22 

, 8 1 Ud 303ii rj 

3 

02 E 0 . 1 . on (030*0 *n* o) .-n 

1 4- A » (VN‘*JN)H 
(N3'N3)H » X 
1 V X - (N3'N3)H 


1 

I 


% 



■61 





310 CONTINUE 
C . 

00 TO 330 

c; 

C COMPLEX PAIR 

320 WR(WA) - X‘+ P 
WR<EN) » X + P 
WI(NA> ^ ZZ 
WHEN) » -ZZ 
330 EN = ENM2 
CiO TO GO 
C 

C ALL ROOTS FOUND. BACKSUBST I TUTE TO FIND 

t' SECTORS OF UPPER TRIANGULAR FORM 

3^0 IF (INUM.EO.n RETURN 
NORM-O.OEO 
K 1 
C 

DO 3G0 I * 1 , W 
C 

DO 3S0 J ^ K. N 
350 NORM‘NORM<^ABS(Hn, JI) 

C 

K - f 

3G0 CONTINUE 
C 

IF (NORM .EQ. O.OEO) GO TO lOOl 
C 

C FOR EN---N STEP -S UNTIL 1 DO -• 

DO 000 NN ^ I , M 
EN - N »■ I - NN 
P -> NR(EN) 

G - W I ( EN ) 

NA - EN -- 1 

IF ( 0 } 7 { <V f GOO < 000 ( 




250. 

2B0. 

251 . 


2B3. 

2G4. 

265. 

2BG. 

2G7. 

268. 


272. 

273. 

274. 
273. 
27(3. 

275. 
270. 
260. 
201 . 
2B2. 


284. 

285. 

286. 
287. 
28B. 
ZSS. 



n n n n on 


GOO 


GIO 

G'/O 

630 


G40 


G50 

700 


M * EN 

H(EN,EW) s l.OEO 

IF (NO .EQ, 0) GO TO 000 

FOR I«EN-1 STEP -1 UNTIL 1 DO -- 
DO 700 II » i, NA 
1 “ £N - II 
W = H (I , I) - p 
R ” H( I .EN) 

IF (M ,GT. MA) GO TO G20 

DO GIO J ^ (1, NA 
R - R + H( I , J) » H( J,EN) 

IF (WI(I) .GE. 0.0E:0) go TO 630 
22 W 
S = R 
no TO 700 
M I 

IF (WKI) .NE. O.OEO) GO TO G 40 
T ^ W 

IF <U .CO. O.OEOJ T - flACHllP • f-JORM 
~ -R / T 
GO lU /nil 

SOLVE REAL EOUATiONS 

X -* 1 1 n , I ♦ n 

Y « H ( 1 1 J , I) 

a * (uiun " p) » (urun - pi <- wki/ * mhi) 
r -I r/i > s " 22 R) / n 

Ma,ewj - T 

IF ( Aesoo ,L£, AOS « 22)) GO TO G50 

HU-H iF.rn - <"R " N It T) / K 
GO TO 700 

H( I H ,LN) ^ (-S Y » r> / 22 . 

CONTINUE 


231 . 

232. 

233. 


293. 

296. 

297. 
203. 

299. 

300. 
301 . 

302. 

303. 

304. 

305. 
30S. 

307. 

308. 

309. 

310. 
311 . 
312. 
3)3. 


315. 

31G. 

317. 

318. 

319. 

321. 

322. 

323. 

324. 





U U u U r. 





GO TO 790 353 

770 M ^ I 35g] 

IP (WKl) .NE. O.OeO) 00 TO 700 3S0. 

23 - Cf'!PLX(-RA,-SA> / CMPDUM.O) 

fi<r.NA) = T3(l) 3G2. 

H(I»£?N) = T3(2) 3G3, 

GO TO 790 3G4. 

C 

t: SOLOE COMPLEX EQUATIONS 

7L’0 X = H(I,I + 1) 3GG. 

Y ^ H( I + l , I ) 367, 

OR = (WR{ I) - P> * (WR( 1) - P) + HI < I ) « WI< I ) - Q # 0 3GB. 

01 = (WR(I) - P) * 2.0E0 » Q 3G9. 

IF (OR .EG. O.OEO .AND. 01 .EG, O.OEO) OR = MACHEP * NORM 370. 
X * ( ABS(N) + ABS(Q) *■ ABSOO 4 A8S(Y) + ABS(2Z>) 

23 ^ CMPLX<X»R-2Z»RA+Q«SA.X«S-Z2«SA-Q«RA) / CMPLX(OR.OI) 372. 
H(I.NA) = T3(l) 373. 

H(I,EN) » T3(2) 374. 

IF ( ABS(XJ .LE. ABG(Z2/ + ABS(Q)) GO TO 785 

H(l4l,NA) (“RA - W 4 Hd.NA) * 0 * Hd.EN)) / X 37G. 

H(I+1/EN) = ("SA - W » H(I.EN) - Q » H(I.NA)) / X 377. 

GO TO 790 37B. 

785 23 = CMPLX(-R-Y*H{ I ,NA) , -9- V4H( I ,EN> ) / CMPLX<ZZ»0) 

ll(I + l,NA> - rj<l) 300. 

H(I + 1,EN) T3(2) 381. 

790 CONTINUE 382. 

END COMPLEX OECTOR 

000 CONTINUE 3B4. 

END BACK SUBGTITUTION. 

OECTORS OF ISOLATED ROOTS 

DO 840 I = 1, N 3B7. 

IF (I .GE. LOW .AND. ) .LE. IGM) GO TO 840 388. 

C , 389. 

DO 820 J - I , M 390. 






620 

840 CONTINUE 

MULTIPLY BY TRANSFORMATION MATRIX TO GIOE 
‘UECTORS of ORIGINAL FULL MATRIX. 
FOR J=N STEP "1 UNTIL LOW DO -- 
DO 880 JJ = LOW. N 
J = N + LOW - JJ 
M = MINO< J . IGH) 

DO 880 I - LOW. IGH 
2.7. = O.OEO 

DO 8G0 K s LOW. M 

SGO ZZ - Z2 <• Z( I .K) HIK.Ji 

Z< I . J ) 22 

880 CONTINUE 


GO TO 1001 

SET ERROR NO CONUERGENCE TO AN 
S EIliLW'.'ALUE AFTER 30 Itl; RAT IONS 

1000 lERR = EN 

1001 RETURN 
C 

C LAST CARD OF HQR2 

END 

END OF FILE 



391 . 

392. 

393. 


335. 

398. 

397. 

398. 

399. 

400. 
401 . 

402. 

403. 

404 . 

405. 
4 0 S . 

407. 

408. 

409. 

410. 


413. 
414 . 


418 


SAMPLE INPUT FILE FOR “SAF 




37.0 


05 

05 


270. 


0 . 

0 . 

0 . 


0 . 
0 . 


HAST ARH-1 OLD DATA SYMETRIC ROOT LOCUS 
i it 1 2 5 .25 o 

15000. 

15 

- . 2000000000000E + 01 0. O. 

‘ 0. 0. 

0 . 0 . I.) . 

0 . 0 . 

9' ' . 1 OOOOOOOOOOOOE + 0 I 0. 

0. I.t. u, 

0. 0 . o. 0 . 

*•* • • 0 . 

. lOOOOOOOOOOOOE+01 -.87202090000001-^05 - . 590GOOOOOOOOOE+03 0. 

O . 0 , , 

*•* • 0 . 0 . 

• i;> . 0 . 

*•* • • 0 . 0 . 

. lOOOOOOOOOOOOE+01 0. 0. 0. 

0 . 0 . 0 . 0 . 

0 . 0 . 0 . 

.2250000000000E + 04 •, 

.2400000000000E+03 0. n 

, • • V • 

'>• 0, 

. • > . ,1 ^ 

V • '•* • • 0 . 

. lOOOOOOOOOOOOE+Ol 0. 

0. 0. 

. i;t . 

.2250000000000E + 04 -.297G191000000E+05 

•. 1G32200000000E+03 - . 22500000000o0E+07 - . 3000000000000E+04 0. 

0. 0. 0. 

*•' • 0 . 0. 

'■> • 0 . 0 . 0. 

- '.v I ', . 

. lOOOOOOOOOOOOE+O 1 0. 0. , 0, 

0 . 0 . 0 . 


1 163G40000000E+0G 


0. 

0 . 


0 . 



<» 


0 . 

- . 1B32200000000E + OG 

- . lOOOOOOOOOOOOE + 03 

0 . 

0 . 

0 . I 

0 . 

o. 

0 . 

12045G3G00000E^<■/6 

0 , 

0 . 

0 . 

0 , 


0 . 

-.21 5G3S4000000H+ 1 0 

0 . 

0 . 

0 , 

0 . 

0 . 

0 . 

0 . 

1 59133SG32000E-* J 0 
-.12544000 0 0 0 0 0 F + < '< 5 

0 . 

0 . 

0 . 


. 2250000000000E+07 - . 297G 1 9 1 OOOOOOE+OG 
- . 2880000000000E+07 - . 33G4000000000E+04 

0 . 0 . 

0 . 

0 . 0 . 

• 0 , 

. lOOOOOOOOOOOOE-^01 o. 


<■' . 0 . 

. I OOOOOOOOOOOOE+0 1 0 . 


o , 


1GG0500000000E+07 - .219S428S58000E+0a 
.2125440000000E+07 . 1 23762GOOOOOOE+04 
IOOOOOOOOOOOOE+03 0. 


0 . 

0 . 

0 . 


0. 

0. 

. I6G05C0000000E+07 

■ ___ 

' . 1 2045G3600000E*OG 

• 1 59 1 33G632000C'* 1 0 

-.21 254 40000000E4 07 


0, 

. 1 SGGOOOOOOOOOtr *05 

0. 


- . 7G7B000000000E+02 

i:> . 

0 . 


(.'1 . 

0. 

I'j. 


0. 

0 . 

0 . 


0 . 

I'j 






0. 

(.) « 

0 . 

. i mOOOOOOOOOOOE+01 


0. 

0 . 
0. 

0 . 
0 . 

0. 

0 . 

0. 

0. 

• 

0 . 

. 1 93 4 2 4 4 OOMOriE"^ OG 

- . 5894OOCO00000E+03 


.100 0 0 0 0 0 0 0 0 0 0 E + <:» 1 

0 . 

0. 

0. 

0. 

0. 

0. 

0 . 

0. 

Cl . 


1*1 

0. 

0. 

c> . 

• 

<1) . 

0 . 

0. 

0. 

0. 

0. 

0, 

0. 

0 . 

0 • . 

0 • 

f'f 

0 . 

. 1 934000000000E+OG 

0 . 



219G42S95a000E+08 

1237G2GOOOOOOE+04 


. 8720209000000E'< 05 


c: . ' 

C * 


*^Tl »>• 

< U) 


OF POOR 



/ 



i 


.825 

T 

•D2039179144333E-01-. 11654315425123E-01-. 125S2a343QSa26E+00 . 12991GlG:531590E+00 
^ 1G731513315278E+00-.549065504G791SE-01 . 1 1 2322421 90525E+00-. 33951509821 190E-01 

.25GG125039S087E+00 . 2372Q743375749E+00 . 101 iai6734603G£-02 

. 58G2 1 245G390G0E+00 - .2511 0777G47203E+0 1 . 1 2924344ea0389E+01 - . 1 203 GGG 3230 G 93 E +01 
.H35S97320041G2F*-«00 . 1 6822540652744E + 0 1 - . 1 703G10ia944G3E+01 . 1 41B9123542322E+01 
- .9284G729342097E+00 . 1 6730890244924E+00 . 70324919GS047GE+0H 

-.201 02552364 770E+01 . 1 347980 1 4 1 5545E+0 1 - .28893S200a7497E+0 1 ~. 7563 1 83813568 lBE +00 
.5355157287G444E-02-. 1021 170Q5744I lE+01 . 1 1 01 83B040573GE +0 1 - . 244 1 1 1 7 1 25750GE+0 1 
-. 45899a23654936E^0I . 358740 1 3 1 77045E+00~ . 1 1552424058953E+01 
.328482653 5305aE^ 01 . 30501 324024G54E + 00- . 394 1 60946G45QGE+00- . 53470 1 66785280E+00 
. 7343862965058BE+00 . 95398101 539840E+00 . 4070a346G37868E-0 1 - . 72 1 B27827B2854E-0 1 

I . 15eeG02189a916E+01-. 1G1G4745279093E+00 . 534701 S94872 5 lE +00 

I - .246833320 15949E+01-.G063S874363G25E+00 . 559451 27470874E+00 . 1 91 9 1 6329 1 9 1 52E+00 

J • .527780837085346 + 00-. 576649035429 1 4E+00--.52223372D94092E-01 . 27506987432007E + 00 

I % -. I695S8648B9223E+01 . 1 534 758909694GE+ 00- . 29 1 97G2B878675E+00 

I -.55199340125432E+00-.5023S2110G1420E+00 . 71 023294 70 1 1 55E+00 . 477084 3208G3GOE+00 

I I -.4795093017G7G4E + 00 . 5540 1 I 11 1 5 1 1 33E + 00- . 1 57092358 1 0732E + 00 . 75 1 B5332B4aa83E-0 1 

■ “ .S517393G262I30E + 00-.92a9I3GOg4'«570E-01 . 8206 1 74G 1 83354E+00 

“.9849a70439G80GE»00-, 670968596559026-01 . 2307929593771 7E + 00 . 8 1 B5 1 303 J 50363E+00 
■}' - . 90595 1469B6539E+ 00 . 1 2797302 1 936 1 06 +00- . 25009025732607E+00 . 5 1 42426975GB72E+00 

-.99314303749173E+00 .262855763099886+00 . 24GB0S3G427240E+00 

- , 43482670248692E+U0 .540667812559816+00- .4424 1 0 1 395 t 027E+00- . 384 355997926G9E+00 
. 142282106532;j36+'.o . 5 t I 9:!600855;?4 1 F +>.,, . 1 (/.) J 606983 1 443E+0 1 -. 27097624 36 1 9G2E + 0 1 
3 12381 03462393E + ./ 1 ", 38694 7830485556 + 00- .261 361 45461 074E + 00 
-.730500910317876+00 ,115053151552336-01 .669185037430356+00 . 24 1 560 1 2804542E+00 
- .342e346B442782E+00 .7667051 1740508E+00 . 88394703731092E+00- . 1 7048 1 7G80B006E+b 1 
-.238757420434486 + 01-. 40l43a077ZlG05E+00 , 5481 98728352 1 9E+00 
i . 258B7043032005E+ 00 - .29191 967602 1 736+00- . 38 19090091 7G04E+00- . 80921 1 29390 1 1 OE+00 

I . 7345B594624844E+00 . 1 5 1 4 4 63547545 1 E+00- . I 0467B2BB 1 7440E+0 1 -. 559347793 1 579 1 E+00 

! -.S3569905G29GG4E+00 . 37797 1 52 1 626936+00- . 60 1 39778707298E+00 

-.354342424902826+00- , 123363294726746 -01-.37303G58755999E+00 . 32397G09793043E+00 
- .547544577916826+00-, I 779040205SB07E+00 . 6350774 1 G36499E+00- . 1 0993986800G09E+0 1 
-.520082039475556+00 . 2929362 1 564 1 40E+ 00 ,322288515854336+00 
-. 2831 62544359036+ 02 -.3692225320274 46+0 l-,3G55403S997203E+01 . 


2 08 6 3 2 2 7 3 1 300 2E + 0 1 










.23410501323667E+01 .34748537215446E+01 .50723BB4273585E+01-. t094306GSBG8lBE+01 
. lS90233B09743aE+01-.4Q312Q22837S27E401-. 1G37473B44IB63E+01 
•.87174G7303379BE+03 . 991 4789421 1 E04E+0 1 . 1922084934015BE+02-.4BG00103690884E+02 
.5255529537528SE+02 .S5457822999557E+01- . l0S9USa275147E+02 . 18073525856157E+01 
.47459335700241E+01-. l0BGieB45312B2E+02 . SeGSe75S88564SE+00 

. 4410S495890S9eEff02-. 734828189701 1 1E+01-. 170798S434034BE+02 . 484853aO20Q144E+O2 
4 180aB4032G354E+02-.43472723397928E+01-. 29025708391 103E+02 .21G157429B04S7E+01 
.535B294075ai84E-01-.5917G413959B17E+02 . HB745045B1 756E+02 
-.2782475792791 lE+03 . 1 BG770a9843947£-* 02 . 31 509S98S14B0SE+02- . 1 144 1377471 712E+03 
. 11 1GQ34460B313E+03 . 4018G091 2838 1 BE+0 i 242806969591 21E+02 . 691 74521 658 123E-01 
.35G007312G141 1E + 0:1 ,3G3013B0315439E + 02- . 1 71 1 2351 187734E + 02 

.740555l31B2B24E+02-. 1 379564 1 1 43 IG7E + 02- . 2502051 227 137SE+02 . 99 1 6803767 1 44 1 £•♦ 02 
- .96440927084700E+02 . 793315S0042085E+00 .2B34a3o3843109E+02-.S722302B0S0397E+00 
.464B7135474532E + 01-.3287305141830BE»-02 . 1 454313773357ZE + 02 

.66aB9540a99283E+02-.5B209442405572E+01-.G2799B7130219lE+01 . 1305G769G92B40E+01 
-.214BB28a422260E+01 103B9941389754E+02-. 1344S5273374B0E+02 . 2191090G5S8037E+01 
1048Gaa5851200E+02 . 274G7G13GB413GE+02- . 815G3492557343E+01 
.39322G182B4221E+03-. 147531879030359+02-. 140446 1925S79SE+02 . 1 35887 1 7047707E+02 


•••. 12875303754B39E + 02-. 109932415a6G04E + 02-.4a9604195G3l77E+02 .30B77047842612E+01 
-,a932£974537793e+01 . 1 70B069304 1 7 1 5E+02 . 159G0474B37GS9E+01 

~. 4 2O32OG3G08222E + O2 . 529 1 91 427774 1 2E-t 0 1 ~ . 38323305457B88E+0 1 . 1 72G9297355304E+01 
-. 659679186 11B50E+00-. 645 1BG89BB3G7GE+00- .251 7743057289BE+00-. 431 057G7851042E+01 
-.690B4a3B2a5713E+01-.8SG999342B6578E+01- .a34785252497BlE+01 
. 1 1B7B4575574 20I: mi:; . 5n338735)n5n4B2E m... - . (;ii>i;j42e7 1 B7564E + 0 1 . 424 7 1 838 1 B2B85E+ 01 


- .S25527G2873968E-M.ll”. 107371 70366 134E m:i 

•• . 1G833334417793E+02 . 3G57G55 1 592823E + 02 

- . 402277483610 18E + 03 . G50334B3503G75E-t0 i 
.37GS264225GS66E+01 . 15G1G14362334GE+02 
. 2 137SG49759334E + <-2 990845586081 46E + 02 
. 1 304328535369 1E+03-. 10300579229 1 45E+02 

■ . 187332 1 7323978E + 02-. 1 B5050337G0058E-M,t2 
- . 131534G5477G71E+02 . G0O773 1 3 1 0780 1 E+0 1 
- . 5484 4296202973E-* 02- . 4 1 4 789 1 75B 174 1.E-* u 1 


14 193SS70672S4E+02-. 39794 104653509E+00 
14G2aB330183G0E+02 

.G35758445B027aE+01 . 70327528G98101E+00 
.3002435G332322E+02-. 10222G73905792E+02 
. 1 171B3140B41 19E+02 

• . 1294I7a47S3901E+02 .2033441 1014025E+02 
- . 25069G37 1 69307E+02 .581 G3G60Q5828 1 E+00 
-.5573900922G153E+01 

- .64323501 424803E+01 . 2706350359a36BE+02 


-.2908571 1 24272 1E+ 02 - .86497331 64 4 IS/E” I. ! 
.2250G703373G98E-MI I - . 349066 1 B8653a0i;: mI I 

391379387737736 M,..( • . l 1 7 1 93863 1 m,;: 


24 11 377267700GE-* 02 •• . 256370 i 87ea490E-* 0 1 
2390409357fJ243i: + 01 

2.?4 2738'l 9600616 + 02 . 1 74883 1 7932SG9E+ 03 


Si 

a 



. 19147263G72297E + 03 
. 12943504932573E + 03 
. 144210lBaB7905E^03 
.848 17^1 466883 14E + 02- 
. 1304G1527ZG081E+03 
. 13G04205427455E + 04- 
.617805G9634494£i+02 
. 740B70B64GG873E+02- 
.5101B275664745E+03 
.8332G4392494y8E-*-02- 
.839559444026286+02 
. 27305530225092E ■T 08- 
. 172G4375578391E+02- 
. 2 1 6GG940007385E+ 0 1 - 
. 1 5523a6GlG655eE+04- 
. 8796820441 4007E+ 0 1- 
.aG0a87G071475GE+Ol 
. 19224108495903E + 03 
. 15775115193230E+02- 
.215591379262176402 
. 515G4536G94907E+03 
.32601 i573G7370E+02- 
.3G7777097303^lE^01- 
. 17476351527783E+04- 
. 799094 7G8G95G9E i 02 
. 35420782874 13:3E^|.l2 
.59061 879 1258G3E+03- 
. 3'^G&4G52O34G70E+O2- 
. 4B990052B435GGE+02 
.205B425440aG71E + 02 
.31 161537477131E+02- 
.42442003207834E+01 - 
.8719800B824936E+03 
. 4939503479432 lE+02- 
. 4 95209500 1 0003E + 0 1 - 
. 731264 1 45:,'2003H.4 02- 


. 1938300232G155E+01 
. 62G352532085O3E+O3- 
. 5077005 1431 440E+01- 
. 1 87162S6087002E+02 
, 563 12342897 197E+ 03- 
.285B59162B5834E+01 
. 17575822G73003E+02 
. 299249252983GGE+03 
. 1003815Q30353gE + <U- 
. 14755G918B5287E + 02 
.3G02567391G625E*03- 
. 335840560 13845E+01- 
.2B5394740G75B5E+01 
,3191 882 ' 899603E + 0 1 ~ 
. 900 1 B805027577E+00- 
.G7053BB57GG07BE+0I- 
. 235492B209B8B2E+02 
.4552850218751BE+01 
.5S105759O37244E+O1- 
.78157595830413E+02- 
.3297201Bai9310EM;U 
.91 44035048032 lE+01- 
. 17332582422323E + U2- 
.374705S2927037FO01 • 
,5 i 1 5337028 1333E fOi 
. 1 75374 ;9305902F: + 03- 
. 145Ga247294605E+01- 
. 13944588975457E+02- 
. 167734213624502E + 03- 
. 44742053770105E+01 
.7B13319G6B3B18E+00- 
. 14122121370427E+01- 
. 37739309233401 E . ’ 

.275265290997 loft '• 
.6 164721 B4744G7E+01 
. 1 08!369023g6G30E I 03- 


. 11B3S1B0374785E+O3 
. 14220BB9552270E+03 
. 70G0977055 1 4G4E+0 1 
. 1 1376O28087B7OE+O1 
. 104B461B271959E+03 
.808375963 12301 E+00- 
. 1814BGG0941254E+01- 
. 50581 039434504 E+02 
. 42962042481220E + 01 
. 242325 18206760E + 02 
.67134500147037E+02 
.221BG2B029154BE+01 
. 46401 7S5BB16B9E+01 
. 6B2G99495993B5E+01 
.2B978B81G48732E+01 
.289804430034B5E+02 
. 1044710B389754E+02 
.gG978B1947B809E+00- 
. 1 l6B6eS7470484£ + 02 
. 2622S09G0B4333E+02 
.2B2B911051B515E+O1- 
.340851 55 1B3512E+02 
.94202B33464343E+01 
-.9514B205444707E+01 
. i!3' < 1 0153B84B0aE-i 02- 
•.434915B84S5755E+02 

- .3G537G30156258E + 01 
-.215004G31501B2E+02 
- . 50971 753777993E+02 

. i44499308G9777E+01- 
•.43053777830211E+01- 

- .87880 177026054 E+00 
. 1060409803301 1E+ 01 
.33165G52037777E+01- 

• , 2 1 9524Z3424222E-1 0 1 
■ .20 1 98454 158909E + 02 
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-78. G4 

0.00 

73.64 

-33.60 

60.06 

73.90 

-33.60 

-68.06 

75 . 90 


PREOUENCy 
CYCLE /BEC 

DAMPINO FACTOR 


1733.67 

l.bOOO 


587.47 

-.0708 


597.47 

.0708 


422.47 

1 .0000 


l83. 71 

1 . 0000 


26 1.68 

f.OOOO 


201.82 

!344 


201 .82 

.1944 


209.58 

l.OOOO 


102. 03 

.3279 


1 02 . 08 

.9279 


105.08 

.0031 


105.00 

.0031 


86.92 

.7529 


66.92 

. 7528 

O o 

88.28 

-.2471 

■VI 38 

88.20 

-.2471 

poo 

73.65 

.0375 

70.65 

.0375 

73 . 74 

.0086 

33 f- 

73.74 

. 0-386 

ny 

66.44 

66.44 

,0085 

.0095 

6<>.34 

.2943 

{■: 

60.34 

.2948 


49.35 

.0416 


49.35 

.0416 


40.76 

.4875 


40.76 

.4875 


40.05 

.0057 


40.85 

.0057 


32.31 

.0148 


32.51 

.0148 


22.68 

1 . 0000 


23.68 

.1320 


23.68 

.1320 


20.23 

.0751 


20.. 13 

.0751 


16.26 

.0138 


16.26 

.0190 


12,61 

1.0000 


12.17 

.4427 


12.17 

.4427 




a 



-33. Q1 

3^.03 

52.23 

-39. et 

-34.05 

52.23 

-2.00 

0.00 

2.00 

-3s.3a 

292.73 

295. 30 

-30.33 

-232.79 

295 . JO 


THE GAIN FACTOf? - ,500 

i 


:iGENVAl.ueS 

COMPUTED. ERROR CODE* 0 


REAL PART 

IMAGINARY PART 

FREGUENCY 

RAD/SEC 

RAD/BEC 

RAD/SEC 

-J08I2.44 

0.00 

10812.44 

203.73 

3718.91 

3726. 2Q 

283.70 

-3716.01 

3726.26 

-2634.35 

0.00 

2634.35 

il4S.30 

0.00 

1145.90 

-1870.04 

0.00 

1670.04 

-230.51 

1230.03 

1281.58 

-230.51 

•1230.03 

1261.58 

-1 176. OS 

0.00 

1178.06 

-880.79 

245.67 

731 .29 

-830.79 

-245.67 

731.29 

- r , 02 

655.30 

053.30 

-2.02 

-655.30 

655 . 30 

135.94 

533.71 

550.75 

135.94 

■533.7! 

550. 75 

-413.30 

364.61 

555.65 

-419.30 




•Ulu .G‘J 

40/. ‘/.i 

-10.54 

•496.69 

-iO/.MO 

-3.05 

460.30 

460. J7 

-3 . '■'5 

-4 GO, 36 

46'.. 3 / 

-3.91 

414.34 

4 1 4 . JC 

-3.91 

• 414.34 

414.30 

-08.94 

370.00 

.181 .40 

-99.94 

**370.88 

:joi . 40 

-12.65 

JO 7. 24 

JO 7. 51 

-I2.S5 

"307.24 

30 7.51 

-1.44 

254.77 

254 . 70 

-1.44 

-234,77 

254.70 

-87.70 

227.04 

243.98 

-07.70 

•227.04 

243.90 

-3.71 

202.77 

202 . 00 

-3.71 

-202.77 

202. 3«* 


A 


t 



a. 37 

.7504 

0.37 

.7564 

.32 

1.0000 

47.35 

.1300 

47.35 

.1300 


FREQUENCY 

CYCLE/BEC 

DAMPING FACTOR 


1733.07 

1 . 0000 


507.47 

-.0706 


507,47 

-.0700 


422.39 

1.0000 


183.73 

-1.0000 

o o 

289. OS 

l.OOOO 

t; » 

202.28 

. 1091 

■0 ?3 

202.20 

.1091 

O 

190.57 

1 .0000 

o 

117.25 

.9413 

?3 {■• 

117.25 

.9419 

••o 

Cl ‘••I. 

105.07 

.0031 

105.07 

.0031 

1 1 

80.31 

-.2 109 

i..' 

08.31 

-.2489 


03 . 09 

.7548 

-L L i 

85.09 

.7546 


73.60 

.03 73 


70.69 

.0373 


73.02 

. 0*708 


73.02 

.0008 


66.44 

,0094 


06.44 

.0094 


61.15 

.2332 


61 .15 

.2332 


49.31 

.0411 


49. ,11 

.0411 


40.85 

.*7*757 


40.85 

' .0*757 


39.12 

.3598 


33.12 

.3590 


32.52 

.0183 


32.52 

.0183 






-138. 1 1 

0 , 00 

1 3U . 1 I 

22.15 

1.0000 

-23,01 

151.30 

153.12 

24.35 

.1503 

-23.01 

-1S1.3B 

153.12 

24.55 

. 1503 

-5.53 

135.30 

133.41 

21.71 

. 04 1 0 

-5.55 

-135.30 

;35.41 

21.71 

.0410 

-1.34 

101.42 

101.44 

16.26 

.0192 

-1,34 

-101.42 

101.44 

1G.26 

.0192 

-BO. GO 

0.00 

60. GO 

12.92 

1.0000 

-37. B4 

‘ 54. G2 

GG.45 

10. G5 

.5695 

-37.04 

-34.62 

66.45 

10.65 

.5633 

-23.72 

36.35 

4B.52 

7.76 

.6126 

-23.72 

-33.35 

48.52 

7.76 

.6126 

-2.00 

0.00 

2.00 

.3? 

1.0000 

-33.23 

292.73 

235.30 

47.35 

.1300 

-33.33 

•292.79 

235.30 

47.33 

. 1300 


THE GAIN FACTOR « 1,000 

o EIGENVALUES COMPUTED, ERROR CODE* 0 


REAL PART 
RAD/SEC 

IMAGINARY PART 
RAD/BEC 

FREQUENCY 

RAD/SEC 

FREQUENCY 

CYCLE/6EC 

DAMP I NO FACTOR! 

10612.44 

0.00 

10612.44 

1733.67 

1.0000 

2S3.B0 

3716.90 

2726.25 

537.47 

-.0706 

263.60 

-3716.90 

3726.25 

567.47 

-.0708 

-2633.39 

0,00 

2633.39 

422.24 

1.0000 

-1737.33 

0.00 

1737.35 

278.57 

1,0000 

-227.27 

1247.02 

1267.56 

203,24 

.1793 

-227.27 

1247.02 

126 V.r.U 

■ 203.24 

.1793 

1146.22 

v.u'. 

1 1 4 u . 2 2 

1U3.76 

-1.0000 

-915.25 

337.20 

975.39 

156.33 

.9363 

-913.25 

-337,20 

9 75.39 

156.39 

.9383 

-603.33 

0.00 

603 . -39 

120.01 

1.0000 

-2.03 

655 . 20 

C55.20 

105.06 

. 003 1 

-2.03 

-655.20 

655.20 

103.06 

.0031 

-427.67 

367.26 

56.1 72 

90.39 

. 7586 

-427.67 

-3B7.26 

5I5.J.72 

90,39 

.7366 

135.74 

534.14 

551.1 2 

63.37 

-.2463 

133.74 

-534,14 

551.12 

06.37 

. -.2463 . 

-18.12 

,497.05 

497.38 

79 . 75 

.0364 

•16.12 

49 7.03 

49/. 30 

79. 75 

.0364 

-.75 

4G0.G1 

■160. G1 

73.65 

.0016 

-.75 

-460.61 

460.61 

73.65 

.0016 

•3,84 

414.33 

414.34' 

06,44 

.0093 

-3.34 

•414.33 

114,34 

• 66.44 

.Oii93 
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-57.35 

373.95 

378.32 

-57.35 

-373.95 

370.32 

-12.03 

306.48 

306.72 

■'12.03 

-306.48 

306.72 

-1.43 

254.76 

254.76 

-1.43 

-254.76 

234. 7G 

-54.07 

237.23 

243.31 

-54.07 

-237.23 

243.31 

-5.48 

‘ 201.22 

201 .30 

•5.48 

-201.22 

201.311 

-136.03 

0 . 00 

136.03 

-28.75 

161.53 

164.07 

-28.75 

- 181.53 

164.07 

-5.78 

142.03 

142.15 

-5.78 

-142.03 

142.15 

-1 .91 

lOl .43 

101 .45 

-1.81 

-101.43 

101.45 

-02.01 

0.00 

82.01 

-45.10 

51.63 

68.61 

-45.18 

-51 .63 

68.61 

-16.33 

33.13 

37.21 

-16.33 

-33.13 

37.21 

-2.00 

0.00 

2.00 

-38.33 

292.79 

293.30 

-38.39 

-232.79 

295.30 

THE 

GAIN FACTOR •* 2.000 



lOENOALUES 

computed, ERROR CODE* 0 


REAL Pari 

IH.iti IWi'.KV (>AI»r 

FREQUENCY 

RAD/GEC 

RAD/SEC 

RAD/SEC 

-10812.44 

0.00 

10812.44 

263.81 

3716.89 

3726.24 

263.81 

-3716.89 

3726.24 

-2631.41 

0.00 

2631 .41 

-1812.64 

0 . 00 

1812,64 

1146.86 

0 . 00 

1146.86 

-208.12 

1282. 85 

1278.88 

-208.12 

•1262,85 

1279.88 

-1014.27 

553.18 

1155.32 

-1014.27 

-553.18 

1155.32 

-708.98 

0.00 

709.30 

-2. 10 

655 .06 

655,06 

-2.10 

655.06 

655.06 


60.60 

.1516 

50.66 

.1516 

49.10 

.0394 

43.10 

,0384 

40.85 

,0058 

40.65 

.0056 

33.01 

.2222 

33.01 

. 2222 

.32 . 28 

.Yin 

32 . 28 

.02 72 

21 .81 

1 . ocioo 

26.31 

, 1 752 

26.31 

.1752 

22.79 

.0407 

22.73 

.0407 

16.27 

.0180 

16.27 

.0188 

13.15 

1 .0000 

11.00 

.6565 

1 1 . 00 

.6585 

5.37 

.4549 

5.97 

■ .4549 

.32 

1.0000 

47.35 

.1300 

4 7 . .15 

.1300 


I 



-p •••/ 


FREOUENCV 

CYCLE/BEC 

1733.67 

507.47 

537.47 
421 .82 
280.64 
! 8.3. 88 
203.22 
203.22 
tas.24 
1E3.24 
113. 84 
103.03 

' 105.03 


OAMPINQ FACTOR 


1.0000 
-.0708 
-.0700 
1 .0000 
1 .0000 
- 1 .0000 
. 1 026 
. 1623 
.8778 
.8779 
1.0000 
.0032 
.0032 
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-432.60 

367.60 

56 7,69 

-432.60 

-367.60 

567.69 

133.53 

535.15 

552.06 

135.39 

-535. 15 

552.06 

-17.32 

487.01 

497.31 

-17.32 

-487.01 

487.31 

2.95 

457.65 

457.66 

2.95 

-457.65 

457.66 

-3.75 

‘ 414.09 

414.11 

-3.75 

-414.09 

414.11 

-9. 24 

361.74 

361.66 

-a. 24 

-331.74 

361.06 

-6.17 

305.81 

305.97 

-6.17 

-305.01 

305,87 

-31.78 

275.18 

277,02 

-31.70 

-275.19 

277.02 

-1.38 

254,80 

254.81 

-1.38 

-254.80 

254,81 

-2.68 

193.69 

198.71 

-2.68 

-198.69 

138.71 

-t34.84 

0 , 00 

134.84 

-40.06 

164.96 

169.94 

-40.86 

-134.98 

189.94 

-7.89 

144.33 

144.54 

-7.89 

-144.33 

144.54 

-1,83 

101.44 

101 .48 

-1.83 

-101.44 

101.46 

-82.89 

0.00 

82.89 

-47.76 

51.01 

69.69 

-47.78 

-31.01 

69.89 

-9.70 

26.41 

28.14 

-9.70 

-26.41 

.'H. 14 

-1.93 


1 .33 

-38,33 

232.73 

295.30 

-38.33 

-292.70 

295.30 


16 Z6R0S A8E 

LOCATED at: 

REAL PART 

IMAGINARY PART 

RAS.'SEC 

RAD/SEC 

11034.83 

0.00 

-4793 10 

7047.97 

-4783.10 

-7047.97 

3346.91 

7142.08 

3346.91 

7142.08 

-0401.13 

0.00 





91.02 

.7620 

91 .02 

. 7620 

63.52 

-.2456 

88.52 

-.2456 

73.74 

.0343 

78.74 

.0348 

73.38 

-.0064 

73.39 

• .0064 

66.40 

.0090 

66,40 . 

.0090 

50.02 

.0255 

SB. 02 

.0255 

49.04 

.0202 

49.04 

,0202 

44,42 

,1147 

44.42 

.1147 

40.86 

.0054 

40.86 

.0054 

31.86 

.0133 

31.88 

.0135 

21 .62 

1 .0000 

27.25 

.2405 

27.25 

.2405 

23.18 

.0346 

23.19 

.0546 

16,27 

.0167 

16.27 

,0137 

13.29 

1,0000 

1 1.21 

.8638 

11.21 

.6836 

4.51 

.3447 

4.51 

.3447 

.32 

1.0000 

47.35 

.1300 

47.35 

.1300 


I 







7407.26 

0 . 00 

143.42, 

3539.00 

143.42 

-3330.80 

2950. IG 

0.00 

791.25 

173.16 

791.25 

-173.16 

-66S.29 

0.00 

-2.70 

654.83 

-2.70 

‘ -854.83 

-437.72 

366.90 

-437.72 

-366.30 

141.57 

534. 13 

141.57 

534.13 

... * y , I'j : 

495.71 

-17.01 

-495.71 

-3.4G 

452.61 

-3.46 

-452.81 

-4.08 

414.20 

-4. OB 

414. 20 

-13.26 

303.79 

-13.26 

-309.79 

-1.44 

254.80 

-1.44 

-254.80 

-60.00 

300.06 

-60.00 

-300.06 

-1.21 

200 . GO 

-1 .21 

- 2C»0 . Stj 

-133.52 

0 . 00 

-9.80 

145.05 

-9.00 

-145.05 

-1 .00 

101.45 

-1.00 

101 .il^ 

-03.92 

■'.Of, 

-50.00 

160.33 

-50.00 

-100.39 

-50.00 

50 . 4 1 

-50.00 

-50.41 

.00 

.14 

.00 

.14 

-.16 

0.00 

-36,39 

292.79 

-30,39 

-292.79 
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